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ABSTRACT 

This review addresses the state of research that employs astronomical (remote sensing) observations 
of solids (“dust”) in young circumstellar disks to learn about planet formation. The intention is for 
it to serve as an accessible, introductory, pedagogical resource for junior scientists interested in the 
subject. After some historical background and a basic observational primer, the focus is shifted to the 
three fundamental topics that broadly define the field: (1) demographics - the relationships between 
disk properties and the characteristics of their environments and hosts; (2) structure - the spatial 
distribution of disk material and its associated physical conditions and composition; and (3) evolution 
- the signposts of key changes in disk properties, including the growth and migration of solids and the 
impact of dynamical interactions with young planetary systems. Based on the state of the art results 
in these areas, suggestions are made for potentially fruitful lines of work in the near future. 

Subject headings: circumstellar matter — planetary systems: formation, protoplanetary disks — dust 


1. MOTIVATION 

The focus on origins , our cosmic context, is innate to 
astrophysics. The past two decades have seen remarkable 
progress and profound, renewed interest in the subject, 
primarily due to advances in the detection and character¬ 
ization of planets around other stars (see Howard 2013 
and Fischer et al. 2014 for recent reviews). There has 
been a considerable observational investment in studying 
exoplanet demographics, quantifying their diverse prop¬ 
erties (masses, orbits; e.g., see the reviews by Marcy et al. 
2005; Udry & Santos 2007; Winn & Fabrycky 2014) and 
identifying some fundamental trends (e.g., with metal- 
licity or host mass; see Johnson et al. 2010). The re¬ 
turn on that investment is a better understanding of the 
key processes that govern the formation and evolution of 
planetary systems. The observed properties of the exo¬ 
planet population place important boundary conditions 
on theoretical explanations of those processes. 

The analogous study of the protoplanetary disks or¬ 
biting young stars offers a powerful complement to that 
work on exoplanets. Disk observations offer insights on 
the pivotal ‘initial conditions’ available for making plan¬ 
etary systems, as well as the early co-evolution and in¬ 
teraction of planets and their birth environments. Ulti¬ 
mately, the observed properties of both disks and exo¬ 
planets can be employed to inform, test, and refine mod¬ 
els of the planet formation process (e.g., Benz et al. 2014 
review a prominent approach to this problem). 

The mutual evolution of disks and their planetary sys¬ 
tems is terrifically complicated. The best theoretical 
models include functionally limited physics, and do not 
yet successfully predict the measured properties of ex¬ 
oplanets given the inferred parameters of disks. Much 
of the lingering uncertainty is due to our relative igno¬ 
rance of disk properties. Disks are rich in information 
content for planet formation models, but observationally 
and physically messy. Fortunately, recent technical ad¬ 
vances are facilitating new, improved disk observations 
that can drive rapid development in those models. 
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With that bright future in mind, this review is intended 
as an introductory pedagogical resource for the observa¬ 
tional side of protoplanetary disk research. The targeted 
audience is junior members of the community (graduate 
students and postdoctoral scientists). The content is not 
comprehensive; the focus is primarily on fundamental is¬ 
sues and popular topics related to disk solids (referred 
to interchangeably as “dust” or “particles”). Although 
all the key techniques are represented, there is a bias to¬ 
ward measurements at radio wavelengths motivated by 
the recent commissioning of the revolutionary Atacama 
Large Millimeter/submillimeter Array (ALMA). 

The review is organized as follows. Section 2 offers a 
brief historical overview of disk observations, to provide 
some context for the ongoing work in the field. Sec¬ 
tion 3 is a primer for interpreting observations of disk 
solids. Sections 4, 5, and 6 discuss the state of the art 
constraints on the demographics of the disk population, 
the inferred structures and physical conditions in disks, 
and insights on how those properties change with time, 
respectively. Section 7 summarizes some open questions, 
with an eye toward future observational research. 

2. HISTORICAL BACKGROUND 

There is a long history of theoretical interest in pro¬ 
toplanetary disks, particularly the precursor to the solar 
system (the solar nebula). Although there are earlier 
metaphysical conceptions, the familiar “nebular hypoth¬ 
esis” of a flattened, rotating structure as the origin of 
the observed co-planar planetary orbits was well formu¬ 
lated by Kant (1755) and Laplace (1796) before subse¬ 
quent elaboration shifted to astronomers (e.g., Moulton 
1900, 1905; Chamberlin 1900; see the historical reviews 
by Brush 1978a,b, 1981 or Wood & Morfill 1988). As 
more physically motivated theories of the star formation 
process were developed, such a disk structure was found 
to be a natural consequence of gravitational collapse in 
a molecular cloud core endowed with some angular mo¬ 
mentum (Hoyle 1960; Cameron 1962; Cassen & Moosman 
1981; Cassen & Summers 1983; Terebey et al. 1984). 

Meanwhile, taking the existence of such disks as a 
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given, the theoretical study of their physical characteris¬ 
tics and evolution was split into two camps; one focused 
on planet formation (see the classic reviews by Cameron 
1988 and Lissauer 1993, or the historical summary of a 
seminal period by Brush 1990), and the other emphasized 
the processes of accretion and momentum transport (see 
Pringle 1981 and Papaloizou & Lin 1995). These top¬ 
ics were mature enough to accumulate a vast literature 
by the 1980s, well before many astronomers were com¬ 
pletely satisfied with the observational evidence for such 
disks (let alone extrasolar planetary systems). 

The discovery and early characterization of young 
stars, the T Tauri (or Orion) variables (Joy 1945, 1949; 
Herbig 1962) and their massive Herbig AeBe counter¬ 
parts (Herbig 1960), hinted at a close connection with 
circumstcllar material. Influenced by the star formation 
simulations of Larson (1969,1972), the blue excesses, line 
emission (e.g., Strom et al. 1971, 1972), and enhanced 
infrared continua (Mendoza V. 1966, 1968; Geisel 1970; 
Gillett & Stein 1971) associated with these young stars 
were seen as evidence for the circumstellar “shells” or en¬ 
velopes predicted to be remnants of the cloud collapse. 
There was debate on whether these features were pro¬ 
duced solely by hot (~10 4 K) gas (e.g., Strom et al. 1971, 
1972, 1975; Strom 1972; Rydgren et al. 1976; Warner 
et al. 1977), or if cooler (~10 2 K) dust was an impor¬ 
tant contributor (e.g., Cohen 1973, 1980; Cohen & Kuhi 
1979; Rydgren & Vrba 1981, 1983; Rydgren et al. 1982). 
Ultimately, it was recognized that both components are 
required; an emerging model attributed the hot (blue) 
excess to accretion shocks at the stellar surface (Lynden- 
Bcll & Pringle 1974; Uchida & Shibata 1984; Bertout 
et al. 1988) and the cool (red) excess to dust on ~AU 
scales (Cohen 1983; Cohen et al. 1985; Lada & Wilking 
1984; Rydgren 1984; Rydgren & Zak 1987). 

That realization was part of a pronounced shift in the 
mid-1980s away from “shells” and toward accretion disks. 
Rather than being precipitated by a single measurement 
or study, this conceptual evolution was driven by the 
confluence of several key independent lines of indirect 
evidence (most of these were clearly summarized in the 
seminal review by Shu et al. 1987): the observed large 
(linear) polarizations of T Tauri stars were indicative of 
scattering off of flattened dust structures (e.g., Elsasser 
& Staude 1978; Vrba et al. 1979; Bastien 1982; Bastien 
& Menard 1988), which were just then being resolved in 
the near-infrared (Grasdalen et al. 1984; Beckwith et al. 
1984; Strom et al. 1985); collimated bipolar outflows 
were being linked to an accretion disk origin (e.g., see 
the review by Lada 1985); asymmetric (blueshifted) for¬ 
bidden emission line profiles were interpreted as smaller 
scale winds or outflows that are partially obscured by an 
extended disk (e.g., Edwards et al. 1987; Cabrit et al. 
1990); the outburst behavior in the FU Orionis variables 
was being associated with disk instabilities (Hartmann 
& Kenyon 1985; Lin & Papaloizou 1985), forging im¬ 
portant links to the fields of cataclysmic variables (e.g., 
Kenyon & Webbink 1984; Lin et al. 1985) and viscous 
accretion disks (e.g., Lynden-Bell & Pringle 1974); opti¬ 
cally thin mm/radio flux measurements confirmed that 
spherical dust distributions were inconsistent with the 
observed optical extinctions (Beckwith et al. 1990; see 
also Churchwell et al. 1987); and perhaps most influen¬ 
tially, increased access to infrared measurements, espe- 



Fig. 1.— An optical HST image of a nearly edge-on disk, seen 
in silhouette against the bright background emission of the Orion 
Nebula (Robberto et al. 2013). It was direct images like this that 
ultimately confirmed the theoretical ideas about the ubiquity and 
basic structures of disks around young stars. 
image credit: NASA, ESA, M. Robberto (STScI/ESA), and the 
Hubble Space Telescope Orion Treasury Project Team 

dally from IRAS (e.g., Rucinski 1985; Strom et al. 1988; 
Harris et al. 1988; Wilking et al. 1989), was motivating 
substantial modeling developments to explain the broad¬ 
band spectral energy distributions (SEDs) in the context 
of irradiated dust disks (Adams & Shu 1986; Adams et al. 
1987, 1988; Kenyon & Hartmann 1987). 

In the minds of many researchers in the field, this shift 
in thinking culminated in a spectacular “seeing is be¬ 
lieving” confirmation with resolved Hubble Space Tele¬ 
scope ( HST) observations of disks in silhouette against 
the bright background of the Orion Nebula (O’dell et al. 
1993; O’dell & Wen 1994; McCaughrean & O’dcll 1996); 
a favorite example is shown in Figure 1. Around the 
same time, spatially resolved images of scattered light 
from HST (e.g., Burrows et al. 1996; Stapelfeldt et al. 
1998) and thermal emission from the first generation of 
millimeter-wavelength interferometers (e.g., Sargent & 
Beckwith 1987, 1991; Koerner et al. 1993; Hayashi et al. 
1993; Lay et al. 1994, 1997; Koerner et al. 1995; Dutrey 
et al. 1994,1996; Mundy et al. 1996; Mannings & Sargent 
1997) also played key roles in shaping the field. 

Since then, the past ~20 years have seen a remark¬ 
able proliferation of disk measurements that have nur¬ 
tured the growth of a new, observational branch to the 
study of planet formation. It is not unreasonable to say 
that unresolved photometry and spectroscopy from the 
Spitzer Space Telescope accounts for a substantial share 
of the maturation of this field, in terms of demograph¬ 
ics, dust composition and basic structure, rapid develop¬ 
ment of radiative transfer tools, and hints of evolution 
(e.g., Hillenbrand et al. 2008; Evans et al. 2009; Furlan 
et al. 2009, 2011). But meanwhile, spatially resolved 
data have been playing a steadily increasing role in those 
advances. With ambitious new technologies being com¬ 
missioned, spatially resolved observations are poised to 
facilitate transformational progress. 

3. ESSENTIALS OF DISK OBSERVATIONS 

Before delving into the state of the field and its po¬ 
tential future directions, it will be useful to go through 
a basic primer on disk observations and how they are 
interpreted in the context of disk properties. 






Observations of Solids in Protoplanetary Disks 


3 


Protoplanetary disks are dense, geometrically flattened 
circumstcllar structures composed of a trace population 
of solids (initially ~1% by mass) suspended in a reservoir 
of (primarily) molecular gas. Gravity and angular mo¬ 
mentum conservation ensure that disk densities decrease 
with radial separation from the host star, r (e.g., Terebey 
et al. 1984; Lynden-Bell & Pringle 1974); thermal pres¬ 
sure and turbulence determine how the densities fall off 
with height above the midplane, z (e.g., Whipple 1972; 
Weidenschilling 1977a; Dubrulle et al. 1995). The disk 
temperatures are almost entirely controlled by the pas¬ 
sive stellar irradiation of the solids (Adams et al. 1987; 
Kenyon & Hartmann 1987), which have large broadband 
opacities that dominate the heating and cooling rates 
(Ossenkopf & Henning 1994; Pollack et al. 1994). Irra¬ 
diation heating implies that temperatures decrease with 
r. Starlight intercepted by dust in the disk surface layers 
(atmosphere) is re-radiated out to space and deeper into 
the disk structure; this results in a vertical thermal inver¬ 
sion, where the midplane is cooler than the atmosphere 
(Calvet et al. 1991; Chiang & Goldreich 1997; D’Alessio 
et al. 1998). In this basic structural framework (see Sec¬ 
tion 5 for more details), the solids completely dominate 
the opacity budget. Coupling this with the high sensitiv¬ 
ity of broadband detectors, it makes sense that observa¬ 
tions of the continuum emission generated by solids are 
the most efficient probes of many disk properties. 

That continuum has two primary contributors, from 
thermal (re-processed stellar energy, which dominates) 
and scattered (reflected starlight) radiation. The specific 
intensity at frequency v along a given line of sight s is 
intimately related to the physical conditions, structural 
distribution, and material properties of the disk solids 
through the formal radiative transfer equation, 

dd v — p A/) ds, (1) 

where p is the density, is the (absorption+scattering) 
opacity per gram of the solid material, and S v is the 
source function (the ratio of emissivity to opacity, with 
absorption and scattering contributions). The condensed 
notation of Eq. (1) obscures the tremendous complexity 
of the radiative transfer problem. A complete solution 
requires the energy input from the host star as well as 
the three-dimensional distributions of p , k v (and the as¬ 
sociated directional scattering properties), and S u . The 
major challenge is that S u also depends on these proper¬ 
ties (p, etc.). In practice, this complicated feedback 
is handled numerically, often with Monte Carlo simula¬ 
tions (Lucy 1999; Bjorkman & Wood 2001) that treat the 
propagation of photons into the disk, the corresponding 
heating of the solids, and the escape of photons to the ob¬ 
server (e.g., Wolf et al. 1999; Whitney et al. 2003; Dulle- 
mond & Dominik 2004a; Robitaillc et al. 2006; Pinte 
et al. 2006; Min et al. 2009; Robitaillc 2011). 1 

3.1. Thermal Emission 

Although appreciating the intrinsic complexity of ra¬ 
diative transfer in disks is valuable, it offers little intu- 

1 There are several Monte Carlo radiative transfer codes designed 
for or amenable to studying disks available as open-source software: 
RADMC-3D, HYPERION, MC3D, and MCMAX are commonly-used examples 
(the latter two require contact with the developers for download). 
Other options (e.g., MCF0ST) are available for collaborative work. 



Fig. 2.— The SED of the young star GO Tau and its associated 
disk ( black points; data compiled by Andrews et al. 2013). Models 
of the stellar photosphere (gray) and an illustrative dust disk (blue', 
computed for a simple parametric structure using the RADMC-3D 
code) are also shown for reference, along with the Spitzer mid- 
infrared spectrum (red; Furlan et al. 2009). 

ition for interpreting the data. But with some simplifica¬ 
tions, a more pedagogical model that captures the essen¬ 
tial features of the thermal emission can be constructed. 
Consider a disk where scattering is negligible, the struc¬ 
ture is vertically thin (S v ~ constant along s), and the 
material is in thermodynamic equilibrium = B U (T), 
the Planck function at the local temperature). Eq. (1) 
can then be recast in terms of the optical depth, defined 
so that dr v = p n v ds, and integrated to give 

h = B v {T) (1 — e~ T ") (2) 

(see Ch. 1 of Rybicki & Lightman 1979 for a deriva¬ 
tion). In this toy model, the continuum is simply black- 
body radiation weighted by the absorbing column. Op¬ 
tically thick emission acts like a thermometer (A, ss B u 
when r„ 1) for the r~ 1 surface; optically thin emis¬ 
sion probes the product of temperature, column density 
( N= f pds ), and opacity (A, ss when < 1). 

The SED is a basic tool for interpreting disk proper¬ 
ties 2 ; an example is shown in Figure 2. At ~l/xm, the 
SED is dominated by the host star. Thermal emission 
from the disk starts to outshine the star at >2-5 /im, and 
peaks in the mid- or far-infrared. This emission is opti¬ 
cally thick: its luminosity reflects the mean temperature 
and area of the emitting region, and its slope is related 
to the radial temperature gradient (Adams et al. 1987; 
Kenyon & Hartmann 1987; Sect. 5.1). Broad spectral 
features (e.g., Fig. 2 near 10/mi) provide some insight 
on the dust mineralogy (e.g., Cohen & Witteborn 1985; 
Waelkens et al. 1996; Meeus et al. 2001; Kessler-Silacci 
et al. 2005; Sect. 5.3). The turnover in the far-infrared 
marks the transition to optically thin emission. The lu¬ 
minosity at longer wavelengths scales with the product 
KyRyilAjust, where Md us t is the dust disk mass (Wein- 
traub et al. 1989; Beckwith et al. 1990; Adams et al. 

2 The SED shows how energy, a useful physical quantity that 
characterizes the radiative transfer of starlight (and intrinsic en¬ 
ergy) through the disk material, is distributed over frequency or 
wavelength. It is defined as vF u or A F\, rather than the practical, 
but physically meaningless, measurement of flux density, F v or F\ 
(= f F dfl). In a begrudging acceptance of standard jargon in the 
field, the latter is referred to as ‘flux’ in this review. 
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Fig. 3. — {top): The cumulative contributions to F u at four rep¬ 
resentative wavelengths as a function of r for the SED model in 
Fig. 2. The cumulative distribution of Md is also shown (in black). 
Note that the 4.5 and 24 pm curves are not zero at the inner edge, 
due to the stellar contribution, {bottom): A two-dimensional map 
of the regions that emit 80% of the flux (or enclose 0.8 M^). 


1990). The mm/radio SED slope probes the shape of the 
opacity spectrum Beckwith & Sargent 1991; Man¬ 
nings & Emerson 1994; Sect. 5.3), itself a function of the 
compositions, morphologies, and sizes of the solids (e.g., 
Miyake & Nakagawa 1993; Henning & Stognienko 1996; 
D’Alessio et al. 2001; Draine 2006). 

The nature of thermal emission (oc B u ) means that 
shorter wavelengths preferentially trace warmer material. 
Since higher temperatures are usually found closer to the 
host star, the SED offers a qualitative A 4 r mapping 
(this is slightly complicated when the vertical dimension 
of the disk is considered; see Sect. 5.1). Figure 3 pro¬ 
vides a guide to the characteristic emission regions at 
four representative wavelengths, using the SED model in 
Figure 2. In the near-infrared, the emission is produced 
by both the innermost (~0.1 AU) part of the disk and the 
stellar photosphere (in this case with nearly equal con¬ 
tributions). Since the inner disk densities are large and 
the near-infrared opacities are high, this emission is very 
optically thick; it originates in a layer that is a substan¬ 
tial height above the midplane. Mid-infrared emission 
has similar optical depths, and probes material at large 
heights on ~few AU radial scales. Far-infrared radiation 
is emitted over a larger range of radii, and can become 
(partially) optically thin in the outer disk (~tens of AU). 
The low optical depths of the mm/radio emission offer 
unique access to cool material in the midplane. 

In general, it is not possible to make unique, quanti¬ 
tative inferences about disk structures from SEDs alone. 
The reason is simple and fundamental. Disks are de¬ 
scribed by parameters that vary spatially (densities, tem¬ 


peratures) or mark specific geometries or locations (e.g., 
the viewing angle, inner/outer boundaries): constraints 
on such quantities from unresolved data are intrinsically 
ambiguous (see Thamm et al. 1994 or Chiang & Goldre- 
ich 1999 for practical demonstrations). That is not to 
claim that SEDs are useless diagnostics; indeed, sophis¬ 
ticated SED models have made crucial insights on disk 
properties (e.g., D’Alessio et al. 1998, 1999, 2001, 2006). 
But, the ability to alleviate or break some key structural 
degeneracies requires spatially resolved data. 

That said, resolved observations of disks are not triv¬ 
ially obtained. Most nearby disks are located at distances 
of ~150pc: there, the solar system dimensions subtend 
0"5 on the sky, and AU scales only ~0 , /01. A quick exam¬ 
ination of Figure 3 indicates that the thermal emission 
from a typical disk is confined well within the diffraction 
limits of ground and space-based telescopes. Interfero¬ 
metric measurements are the solution. Infrared interfer¬ 
ometry offers unique insights on the inner disk structure 
at exquisite angular resolution (for more details, see the 
comprehensive reviews on the subject by Millan-Gabet 
et al. 2007 and Dullemond & Monnier 2010); however, 
the emission at these wavelengths traces a minuscule 
fraction of the disk mass and extent (see Fig. 3). In¬ 
stead, it is (sub-)mm/radio interferometry that provides 
the most robust information on disk structure over the 
widest range of spatial scales. Contrary to a common 
misconception, observations at these long wavelengths 
are not solely probing the cool outer regions; with suffi¬ 
cient sensitivity and resolution (i.e., from ALMA), they 
will also provide rare access down to r ^fewAU. 

Interferometers measure the Fourier transform of the 
brightness distribution at a discrete set of spatial scales, 
determined by the spacings between array elements and 
their projected motion onto the sky plane as the Earth 
rotates. A complex “visibility” is measured at each scale; 



deprojected baseline [kA] 

Fig. 4. — The 870 /im continuum visibility profiles (normalized 
real components, azimuthally averaged accounting for the viewing 
geometry projected onto the sky) for the GSS 39 (black), VSSG 1 
(red) and SR 21 (blue) protoplanetary disks (Andrews et al. 2009). 
The corresponding images are shown in the top panels, with the 
synthesized beam dimensions marked in their lower right corners. 














Observations of Solids in Protoplanetary Disks 


5 



Fig. 5. — Coronagraphic images from HST showing the resolved 
optical and near-infrared starlight scattered off the surfaces of the 
disks around HD 163296 (left-, Grady et al. 2000) and TWHya 
(right; Debes et al. 2013), respectively. Some additional exam¬ 
ples in different contexts can be found in Fig. 11 and 19. These 
images are reproduced with permission from the authors and the 
American Astronomical Society; © AAS. 

in a coordinate frame with the disk origin at the center, 
the real components are sensitive to the radial variation 
of the brightness and the imaginary components record 
deviations from symmetry. The data are often displayed 
in condensed format, by azimuthally averaging into one¬ 
dimensional profiles that represent how the emission de¬ 
pends on spatial scale (Lay et al. 1994, 1997; Hughes 
et al. 2007): Figure 4 shows some examples. These pro¬ 
files roughly correspond to the Fourier transform of the 
corresponding radial surface brightness profile. Recalling 
some Fourier transform pairs should clarify the behavior: 
a point source in the image domain has a constant am¬ 
plitude in the visibility domain, a compact source has an 
extended visibility profile (and vice versa), and a ring- 
shaped source (Sect. 6.2) makes oscillations like a Bessel 
function. The visibilities are usually also inverted to con¬ 
struct an image (top panels in Fig. 4), but any analysis 
should usually be performed in the Fourier domain. 

3.2. Scattered Light 

While the thermal radio emission is extraordinarily in¬ 
formative, it is not the only disk probe sensitive to a wide 
range of spatial scales. Scattered light, primarily in the 
optical and near-infrared, provides another key resolved 
diagnostic (see Watson et al. 2007 for a review); Figure 5 
shows some examples of scattered light images. There is 
no pedagogical simplification for the radiative transfer 
equation when scattering is considered, since the source 
function is coupled to the radiation field. 3 Starlight is 
scattered by small dust grains with high albedos, having 
dimensions comparable to the spectral peak of the inci¬ 
dent radiation field (~l/xm). The corresponding large 
scattering opacities mean that the observed emission 
originates in the disk surface layers, and can be detected 
over a large range of disk densities (and therefore ra¬ 
dial scales). Such observations are fundamental probes 
of vertical structure (Sect. 5.1) and can measure basic 
geometric parameters (e.g., inclination). Incorporating 
color and polarization data constrains the compositions, 
sizes, and shapes of the scatterers (Sect. 5.3). 

The clear advantage of scattered light measurements 
is that they routinely achieve exquisite angular resolu- 

3 Si, includes the integral of I u over angle, which requires an 
iterative numerical solution to Eq. (1). This is not the case for 
pure thermal emission, where S v depends only on temperature (as 
was assumed in the approximation of Eq. 2). 


tion (~0 , /l or less), making them well poised to identify 
and characterize ~10 AU-scale substructure in disks (see 
Sect. 6.2). The major challenge is contrast, since the host 
star is typically orders of magnitude brighter at the opti¬ 
mal scattering wavelengths. Coronagraphs or differential 
imaging techniques have been developed to mitigate such 
issues. Much of the current focus is on polarimetric dif¬ 
ferential imaging (e.g., Tamura 2009; Hinkley et al. 2009; 
Quanz et al. 2011; Perrin et al. 2014), which exploits the 
fact that direct starlight is unpolarized, but the scatter¬ 
ing by dust grains can induce a significant polarization 
(Potter et al. 2000; Kuhn et al. 2001; Perrin et al. 2004). 
In practice, measurements of scattered light near the star 
are limited by subtraction residuals or coronagraph di¬ 
mensions to ~0.1-0"2 (the “inner working angle”). At 
larger radii, the data are typically sensitivity-limited be¬ 
cause of the diluted stellar radiation field (there is a 1 /r 2 
depletion of scattering intensity on top of any intrinsic 
fall-off caused by, e.g., a decreasing density of scatterers). 

3.3. Synopsis 

• Solids (dust grains) dominate the disk opacity bud¬ 
get; observations of disks are most sensitive to the 
corresponding bright continuum emission. 

• Thermal emission is observed over a broad wave¬ 
length range; when optically thick, it measures the 
temperature in a surface layer, and when optically 
thin it traces the product of temperature, density, 
and opacity (dust ensemble properties). 

• The morphology of the infrared SED is roughly re¬ 
lated to the radial temperature distribution. 

• Spatially resolved data is mandatory to learn about 
disk structures; mm/radio interferometry measure¬ 
ments uniquely probe material in the midplane over 
a wide range of spatial scales. 

• Scattered light imaging is a complementary probe 
of resolved structure, especially sensitive to the ver¬ 
tical dimension and sub-structure on small scales. 

Additional Reading: the first chapter of Rybicki & Light- 
man (1979); the review of some SED basics by Beck¬ 
with (1999); the summary reviews on disk measurements 
with nnn/radio interferometers by Wilner & Lay (2000) 
and using scattered light images by McCaughrean et al. 
(2000) and Watson et al. (2007); the overview of disk 
observations by Hartmann (2008, his Ch. 8). 

4. DEMOGRAPHICS 

Planet formation research is still data-limited: the 
complexity of theoretical ideas outpaces observational 
constraints (although perhaps not for long). As in many 
similar areas of astrophysics, this means that progress is 
made in two complementary ways: (1) basic demographic 
studies for large surveys, and (2) in-depth analyses of in¬ 
dividual (or small groups of) disks. This section focuses 
on the former. Case studies and detailed investigations 
of small samples are addressed in Sections 5 and 6. 

Simple and practically accessible observational diag¬ 
nostics are required to accumulate a sample large enough 
to quantify the statistical distribution of some property, 
or connections between properties, in a disk population. 



6 


Andrews 


Key examples include the use of infrared photometry 
to describe the basic structural evolution of circumstel- 
lar material (Lada & Wilking 1984; Adams et al. 1987; 
Meeus et al. 2001; Robitaille et al. 2006; Evans et al. 
2009) and the characteristic lifetime of warm dust in 
disks (Haisch et al. 2001; Hernandez et al. 2008; see 
the summary by Mamajek 2009). Simple accretion di¬ 
agnostics, like the luminosity of the blue excess or the 
strength of a bright emission line, are also commonly 
employed to probe evolution (e.g., Hartmann et al. 1998; 
Sicilia-Aguilar et al. 2005) and disk-host relationships 
(e.g., Muzerolle et al. 2005; Natta et al. 2006). 

From the perspective of planet formation, the disk 
mass ( Md ) is a fundamental demographic parameter: in 
any theory, the efficiency of planet formation depends 
critically on the amount of raw building material avail¬ 
able. Moreover, Md can be estimated from a mm/radio 
flux measurement (oc Md us t; Sect. 3.1) given some as¬ 
sumptions for (or ideally measurements of) the dust-to¬ 
gas ratio, dust opacity, and temperature. 4 This section 
addresses observational insights on the principal issues 
that affect disk masses, specifically the intrinsic relation¬ 
ship with host mass (M*; Sect. 4.1) and externally driven 
evolution through interactions with the local (Sect. 4.2) 
and global (Sect. 4.3) environment. A brief comment on 
what has been learned about the age-related evolution of 
the disk mass distribution is also made (Sect. 4.4). 

4.1. Dependence on Host Mass 

Most theoretical studies of planet formation assume 
that Md oc M*. The intuition for this simple scaling re¬ 
lation follows from the star formation paradigm (cf., Shu 
et al. 1987), where a protostar accretes mass through its 
disk, which is fed by an envelope reservoir. The simi¬ 
lar shapes of the observed core (envelope) and star mass 
functions affirm this kind of “gravitational” scaling (e.g., 
Motte et al. 1998; Enoch et al. 2006), but it need not 
apply to the intermediary disks. Even if a disk/host 
mass relationship were imprinted at the star formation 
epoch, it may not persist throughout its (potentially In¬ 
dependent) evolution up to the ages when disks become 
observable and form planets (> lMyr later). 

And yet, the demographics of the exoplanet population 
offer compelling evidence in favor of such a relationship. 
Radial velocity surveys have demonstrated that the fre¬ 
quency of giant planets in compact (< few AU) orbits 
around nearby field stars scales linearly with host mass 
(Johnson et al. 2007, 2010; Bowler et al. 2010). Since 
the giant planet formation efficiency should roughly scale 
with the disk mass (e.g., Pollack et al. 1996; Hubickyj 
et al. 2005), this has been interpreted as strong, indirect 
evidence for a link between Md and M*. 

The signature of that relationship was lacking in early 
mm/radio continuum surveys (e.g., Beckwith et al. 1990; 
Osterloh & Beckwith 1995; Andrews & Williams 2005), 
although the sampling in M* was restricted to around 
1 Mq . When extended to disks orbiting very low-mass 
hosts (~0.1 Mq), significantly fainter emission is mea¬ 
sured (Klein et al. 2003; Scholz et al. 2006; Schaefer et al. 
2009). A complete 1.3 mm continuum census across the 
host mass spectrum in the Taurus region confirms that 

4 For a dust-to-gas mass ratio the total disk mass is = 
Md us t (1 -I- C 1 ), or W/ ~ M, luS | /e for a typical case where £ -C 1. 



Fig. 6. — (top): The 1.3 mm continuum fluxes (oc M,j) from Tau¬ 
rus disks as a function of their host masses (cf., Andrews et al. 
2013). Upper limits on F v are shown as faded markings with left¬ 
pointing arrows. Despite considerable scatter, there is a clear cor¬ 
relation present that is consistent with a linear Md oc M* scaling; 
the gray curve shows such a relationship for the 2 Myr Siess et al. 
(2000) model isochrone and the Andrews et al. assumptions for 
linking F v and Md- (bottom): The cumulative distributions for F v 
in four different host spectral type bins. These distributions are 
truncated by the survey sensitivity threshold at the low-flux end. 

the luminosities increase for hosts with earlier spectral 
types (Andrews et al. 2013), as demonstrated in Fig¬ 
ure 6. For a fixed opacity, Andrews et al. found reason¬ 
able agreement with a linear Md oc M* scaling and a 
typical disk-to-host mass ratio of 0.2-0.6%. That said, 
there is substantial scatter around that relationship, with 
a 0.7dex (factor of 5) dispersion in Md at any given M*. 5 
Mohanty et al. (2013) independently verify this behavior 
with a different (incomplete, but partially overlapping) 
sample, and tentatively suggested that the relationship 
flattens (or perhaps even turns over) at high M*. The 
same linear scaling is also consistent with the small sur¬ 
vey of disk masses in the (older) Upper Sco association 
recently conducted by Carpenter et al. (2014). 

It is tempting to link the observed scaling relations be¬ 
tween Md, M», and giant planet frequency. Assuming 
a direct mapping (for solar metallicity), and associating 
the giant planet frequency of ~0.07 (M*/Mq) (Johnson 
et al. 2010) with the most massive disk progenitors, the 
Taurus disk mass distribution (Andrews et al. 2013) im¬ 
plies that giant planets form when Md > 0.03 M*. That 
mass threshold compares well with constraints in the 
solar system (Weidenschilling 1977b; Hayashi 1981; see 
Sect. 5.2), and, coupled with the rough linearity of the 
scalings, is consistent with predictions of the core accre¬ 
tion theory for giant planet formation (Laughlin et al. 

5 Note that this scatter means that some more complex rela- 
tionships between M^ and M* (rather than a simple power-law 
behavior) cannot be easily ruled out (Andrews et al. 2013). 
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2004; Ida & Lin 2005; Kennedy & Kenyon 2008; Alib- 
ert et al. 2011; but see Kornet et al. 2006). However, 
the inferred disk masses are considerably uncertain in an 
absolute sense (see Sect. 5.2-5.3 for more details): aside 
from the notable scatter (Fig. 6), the opacities and dust- 
to-gas mass ratio are poorly constrained, and therefore 
the conversion from observed emission (F v ) to Md may 
also be biased. If the inferred Md values are systemati¬ 
cally under-estimated (by factors of a few), then the al¬ 
ternative disk instability mode of giant planet formation 
may also explain the linear scalings of disk/host/giant 
exoplanet properties (Boss 2011). 

4.2. Effects of Host Multiplicity 

Many stars are members of binary or higher-order mul¬ 
tiple systems (see the recent reviews by Duchene & Kraus 
2013 or Reipurth et al. 2014). The multiplicity frequency 
for Sun-like stars is ~50% in the field (Abt & Levy 1976; 
Duquennoy & Mayor 1991; Raghavan et al. 2010), and 
may be even higher in nearby young associations (Leinert 
et al. 1993; Ghez et al. 1993; Reipurth & Zinnecker 1993; 
Simon et al. 1995; Kraus et al. 2011). Moreover, the dis¬ 
tribution of orbital separations in multiples peaks near 
100 AU (e.g., Raghavan et al. 2010; Kraus et al. 2011), 
comparable to the typical disk size (see Sect. 5.2). Taken 
together, these properties suggest that the presence of a 
companion could play a decisive, general role in planet 
formation and the evolution of protoplanetary disks. 

Dynamical interactions in multiple systems alter the 
structures of the associated circumstellar material: the¬ 
oretical models indicate that tidal forces truncate indi¬ 
vidual disks at r > 0.2-0.5a, where a is the orbital sepa¬ 
ration, and clear the inner regions of circum&mary disks 
at r < 2-5a (e.g., Lin & Papaloizou 1993; Artymowicz 
& Lubow 1994, 1996). Despite such a disruptive envi¬ 
ronment, the mature counterparts of these young mul¬ 
tiples are known to host planets (e.g., Patience et al. 
2002; Raghavan et al. 2006; Desidera & Barbieri 2007; 
Bonavita & Desidera 2007). Observations of the disks 
in young stellar pairs offer key insights on how planet 
formation is affected in multiple host environments. 

Given the predicted outcomes of these star-disk inter¬ 
actions, the observational focus has been on characteriz¬ 
ing how disk properties vary with a (or its sky-projected 
equivalent). Seminal early work by Jensen et al. (1994, 
1996) indicated that Md is diminished for stellar pairs 
with a < 50-100 AU, in quantitative agreement with the 
expectations of tidal truncation models (see also Osterloh 
& Beckwith 1995; Andrews & Williams 2005; Cieza et al. 
2009) . Recent extensions of those studies with improved 
sensitivity confirm these results (Harris et al. 2012; Ake- 
son & Jensen 2014): Figure 7 illustrates this behavior 
in the Taurus region. The closest pairs have lower disk 
masses (but see below) than those with wide separations 
or their singleton counterparts; the median Md scales up 
by a factor of ~three per decade in a. Similar conclusions 
are drawn from infrared excesses (e.g., Cieza et al. 2009; 
Kraus et al. 2012). 6 The frequency of short-period giant 

6 However, as was highlighted in Sect. 3.1, the high optical 
depths at such wavelengths makes it difficult to link the results 
to mass depletion. Indeed, some studies find examples of infrared 
excess emission regardless of the projected host separation (e.g., 
White Sz Ghez 2001; McCabe et al. 2006; Pascucci et al. 2008). 



0.001 0.01 0.1 1 
F v : 1.3 mm flux density [Jy] 


Fig. 7. — (top): The 1.3 mm flux (oc Md) from disks in stellar 
(spectral type <M5) multiple systems as a function of the sepa¬ 
ration between host pairs (cf., Harris et al. 2012; includes data 
from Akeson & Jensen 2014). Symbols are as in Fig. 6. The flux 
distribution shifts to fainter values when the disks are hosted by 
closer pairs, although there is a bimodal distribution at the closest 
separations where circumbinary disks ( orange) can be especially 
bright (massive), (bottom): The cumulative distributions for F v in 
different pair separation bins. The gray distribution is for singles. 

planets with binary hosts is also suppressed (Wang et al. 
2014), consistent with this inferred disk mass depletion 
present at the planet formation epoch. 

However, the separation-dependent pair demographics 
tell only part of the story. Component-resolved mm ob¬ 
servations indicate that the primary disk almost always 
dominates the Md budget, and that often the disk mass 
ratios are much different than would be expected from 
the standard theory given the stellar masses and separa¬ 
tions (Jensen & Akeson 2003; Patience et al. 2008; Har¬ 
ris et al. 2012; Akeson & Jensen 2014). Moreover, the 
resolved sizes of individual disks in such systems show 
at best marginal consistency with theoretical predictions 
for their tidal truncation radii (Harris et al. 2012). Aside 
from these apparent discrepancies, the pair demographics 
themselves highlight an additional (often unappreciated) 
mystery: unlike any other separation scale, there is a 
clear bimodal Md distribution for the closest (< 10 AU) 
pairs (see Fig. 7), with most remaining undetected but a 
few having notably massive circumbinary disks. 

Some of these observations might best be explained 
in terms of how the process of multiple star formation 
directs mass from the natal envelope to specific compo¬ 
nents in the system (i.e., initial conditions; Bate & Bon- 
nell 1997; Bate 2000; Ochi et al. 2005). But there are also 
indications that the standard tidal interaction models are 
insufficient; modern hydrodynamic simulations suggest 
that asymmetries (like eccentricities or warps) are funda¬ 
mental (e.g., Kley et al. 2008; Paardekooper et al. 2008; 
Marzari et al. 2009). Particularly striking are the sys- 
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Fig. 8.— The 860 pm fluxes (ocM^) from disks in the ONC 
as a function of their projected separation from the massive star 
fdOriC (cf., Mann & Williams 2010; Mann et al. 2014). Upper 
limits (3 cr) are marked with gray arrows. Disk masses drop precip¬ 
itously when their hosts are located in the small region (<0.03 pc) 
where extreme UV emission from the Trapezium stars dominates 
the radiation field, as expected from photoevaporation models. 

terns with disks that are misaligned with respect to each 
other (Jensen et al. 2004; McCabe et al. 2011; Jensen 
& Akeson 2014; Williams et al. 2014) and/or the orbital 
planes of their respective hosts (e.g., Akeson et al. 2007; 
Verrier & Evans 2008; Andrews et al. 2010a). 

4 . 3 . Environmental Impact 

Most disk observations are necessarily focused on the 
nearest star-forming regions, which have low stellar den¬ 
sities and few (if any) massive stars. However, most stars 
form in the significantly different environments of dense 
clusters (Lada & Lada 2003; Porras et al. 2003), where 
much larger stellar populations create opportunities for 
disks to be affected by two key external environmental 
factors. First, high stellar densities increase the proba¬ 
bility for “fly-by” interactions, and thereby the tidal dis¬ 
ruption of disks (e.g., Clarke & Pringle 1993; Korycansky 
& Papaloizou 1995; Larwood 1997; Boffin et al. 1998; 
Kobayashi & Ida 2001). Second, the high-mass tail of 
the M* distribution is populated; disks in the immediate 
vicinity of massive stars will be depleted by photoevap¬ 
oration, due to the locally intense radiation field (e.g., 
Hollenbach et al. 1994; Johnstone et al. 1998). 

There is not yet any clear evidence indicating that tidal 
stripping from close stellar encounters plays a significant 
role in setting disk properties in rich clusters. Indeed, the 
expected interaction rate in the nearest massive cluster 
(the Orion Nebula cluster, or ONC) is quite small (e.g., 
Adams et al. 2006; Proszkow & Adams 2009). However, 
the relevant observational signature - a substantial ra¬ 
dial truncation, and thereby mass depletion (e.g., Bres¬ 
lau et al. 2014) - is not yet readily available in the ONC 
or other, more distant, regions. 

There is a wealth of data that probes the photoe¬ 
vaporation of disk material (around Sun-like hosts) in 
the proximity of massive stars, typically in the form 
of ionized (shocked) shell structures associated with sil¬ 
houette disks (e.g., Churchwell et al. 1987; Bally et al. 
1998a, 2000; Henney & O’Dell 1999; Smith et al. 2003). 
Early attempts to quantify the mass depletion in these 
disks with mm interferometry were difficult (e.g., Mundy 


et al. 1995; Bally et al. 1998b; Eisner & Carpenter 2006), 
but eventually bore fruit at higher frequencies (Williams 
et al. 2005; Eisner et al. 2008). Mann & Williams (2009, 
2010) identified a trend toward smaller disk masses for 
hosts located closer to the OB stars in the center of 
the ONC. Figure 8 shows a recent confirmation of this 
trend, which provides crucial quantitative evidence for 
photoevaporation-driven mass loss (Mann et al. 2014): 
stars located inside the region where energetic ultraviolet 
photons from the 08 star 6 1 Ori C dominate the radiation 
field (< 0.03 pc) host only faint (low-mass) disks. 

4 . 4 . Mass Evolution 

By the time they are observable, disks have already un¬ 
dergone some (perhaps substantial) evolution. The key 
signatures of that evolution will be covered in Section 6, 
but the relevant timescales are best considered from a 
demographic perspective. The standard astronomical 
approach tracks how a given property changes between 
large samples that have different mean ages. ' Applica¬ 
tion of this technique to simple disk tracers, like the in¬ 
frared excess (e.g., Haisch et al. 2001; Hernandez et al. 
2007) or accretion rate indicators (e.g., Sicilia-Aguilar 
et al. 2010), suggest that the typical inner disk (< 1AU) 
survives over a ~5-10Myr timescale (but see Pfalzner 
et al. 2014). But such constraints offer little in terms of 
assessing the available timeframe for planet formation. 
For that, it is critical to know the mass depletion rate, 
or how the Md distribution changes with time. 

The current observational information on Md evolu¬ 
tion is very limited. Only the Taurus region (^2-3 Myr) 
has a mm continuum census complete to a sufficiently 
deep limit (Beckwith et al. 1990; Osterloh & Beckwith 
1995; Andrews & Williams 2005; Andrews et al. 2013). 
Continuum surveys are also available for the NGC 2024 
(lMyr; Eisner & Carpenter 2003; Mann et al. 2015), 
Ophiuchus (1-2 Myr; Andre & Montmeric 1994; Nuern- 
berger et al. 1998; Andrews & Williams 2007a), ONC (1- 
3Myr; e.g., Mann et al. 2014), Lupus (1-3Myr; Nuern- 
berger et al. 1997), Chamaeleon (1-3Myr; Henning & 
Thamm 1994), MBM 12 (1 -3Myr; Hogerheijde et al. 
2002), IC 348 (3 Myr; Carpenter 2002; Lee et al. 2011), cr 
Ori (3-5 Myr; Williams et al. 2013), A Ori (5 Myr; Ansdell 
et al. 2015), and Upper Sco (5-10Myr; Mathews et al. 
2012; Carpenter et al. 2014) clusters, although with in¬ 
homogeneous sizes, completeness levels, and sensitivities. 
Comparisons between these samples are difficult, due to 
the strong selection effects related to the demographic 
trends discussed above; biased, small samples can mimic 
or obscure real evolutionary changes in the Md distribu¬ 
tion (Andrews et al. 2013). That said, the low detection 
rates in the oldest samples that have been probed (a Ori, 
A Ori, and Upper Sco, at ~5 Myr) do indicate an overall 
decrease in the mean disk mass. 

4 . 5 . Synopsis 

• Disk mass is the fundamental aggregate property 
in planet formation models. The factors that in¬ 
fluence Md can be studied demographically with 

7 The ages of young stars are not well determined (see Soderblom 
et al. 2014). The benchmarks common in the literature and prop¬ 
agated here can be interpreted with considerable liberty. 
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large mm/radio continuum photometry surveys, al¬ 
though are subject to uncertainties in the assumed 
dust-to-gas ratio, opacities, and temperatures. 

• Disk masses are related to their stellar host masses; 
a roughly linear Md oc M* scaling seems appropri¬ 
ate, although there is a large dispersion (0.7 dex) in 
Md for any given M*. For standard assumptions, 
the median disk-to-star mass ratio is ^0.2-0.6%. 

• In multiple star systems, close pairs tend to host 
low-mass disks. Most of the circumstellar material 
is usually concentrated around the primary star. 

• Disk masses are depleted (by photoevaporation) in 
the immediate vicinity of massive (OB) stars. 


model for a disk structure is based on some simple phys¬ 
ical arguments for the gas , implicitly assuming that the 
solids follow the same distribution (but see Sect. 6.1 for 
some important caveats). The underlying principle is 
that the disk is in vertical hydrostatic equilibrium, 

dP 

< s > 

where P is the gas pressure and g z is the z-component 
of the gravitational acceleration. In most cases, the ideal 
gas equation of state is appropriate and the stellar host 
dominates the potential. Eq. (3) is then equivalent to 

fllup 

dz 


/zttih GM* z d In T 
kT ( r 2 + z 2 ) 3 / 2 dz 


• The timescales over which the Md distribution 
changes are not yet clear, due to incomplete sam¬ 
ples and selection effects. Preliminary results sug¬ 
gest that there is significant depletion within 5 Myr. 

Additional Reading: the general review by Williams & 
Cieza (2011); the topical reviews of disk dissipation 
timescales by Hillcnbrand (2008) and Mamajek (2009). 

5. STRUCTURE 

Although demographic studies reveal some fundamen¬ 
tal properties of the disk population, their reliance 
on easy-to-observe “compound” diagnostics means that 
they cannot tell the whole story. An important comple¬ 
ment is found in detailed probes of individual disk struc¬ 
tures, which comprise a set of more “elemental” mea¬ 
surements. The single most valuable of these is the spa¬ 
tial distribution of mass. The efficiency of planet forma¬ 
tion is directly linked to local disk densities: there must 
be enough stuff in the right places (for enough time) to 
assemble a planetary system from its progenitor (disk) 
material. Such a density threshold in solids is especially 
important in the formation models of terrestrial plan¬ 
ets (e.g., Raymond et al. 2004; Kenyon & Bromley 2006; 
Kokubo et al. 2006) and giant planet cores (e.g., Mizuno 
1980; Pollack et al. 1996; Hubickyj et al. 2005). 

Inferring the disk structure from its associated obser¬ 
vational tracers is complicated: there are elaborate de¬ 
pendences between the densities, temperatures, and ma¬ 
terial properties. Fortunately, different kinds of data are 
especially sensitive to specific aspects of disk structure 
(Sect. 3). When complementary tracers are considered 
together, they can be used to forge crucial benchmarks 
for models of planet formation and disk evolution. This 
section covers the key aspects of disk structures, with 
emphasis on the practical connections between physical 
parameters and data. It highlights the role of vertical 
structure in regulating the thermal budget (Sect. 5.1), 
the radial distribution of mass (Sect. 5.2), and the char¬ 
acteristic properties of the disk solids (Sect. 5.3). 

5.1. Vertical Structure 

5.1.1. Physical Overview 

The aspect ratio of a rotationally flattened disk is 
small, z/r ~ 0(0.1), since material is gravitationally con¬ 
centrated at the midplane. But even this small vertical 
extent has important consequences. The most common 


with G the gravitational constant, p the mean molecu¬ 
lar weight, toh the mass of a hydrogen atom, and k the 
Boltzmann constant. Eq. (4) is generally solved numeri¬ 
cally, although some common simplifications enable an il¬ 
lustrative analytic solution. For a geometrically thin disk 
(z <C r) with a small temperature gradient ( dT/dz « 0), 
the solution to Eq. (4) is a simple Gaussian distribution, 


with 





c L= ( kT r 3 y /2 
Q \ /ttoh GM* / 


( 5 ) 

( 6 ) 


denoting a characteristic scale height (the ratio of the 
sound speed, c s , to the Keplerian angular speed, f2; cf., 
Shakura & Sunyaev 1973) and E representing a surface 
density (a boundary condition of the integration). 

Simple energy arguments based on their observed bolo- 
metric luminosities indicate that disks are primarily 
heated by stellar irradiation (Adams & Shu 1986; Adams 
et al. 1987). With their large, broadband opacities maxi¬ 
mized near the peak of the stellar spectrum, small (^/un¬ 
sized) grains are the most efficient conduit for that irradi¬ 
ation energy. Starlight absorbed in the disk surface lay¬ 
ers heats those grains, which then re-emit some of that 
energy deeper into the disk interior to warm the mid¬ 
plane (Calvet et al. 1991, 1992; Malbet & Bertout 1991). 
This external deposition of energy produces a thermal 
inversion in the disk atmosphere (T increases with z) 
and modifies the vertical density distribution (cf., Eq. 4; 
Chiang & Goldreich 1997, 1999; D’Alessio et al. 1998, 
1999; Bell et al. 1997; Bell 1999; Dullemond et al. 2001, 
2002; Malbet et al. 2001). All else being equal, the heat¬ 
ing depends on the irradiated area of the surface layer. 
If H(r)/r is increasing, the tilt of this layer toward the 
star increases that area: such flared disks intercept more 
starlight than their flat counterparts, and therefore have 
warmer temperatures (at a given r; Kenyon & Hartmann 
1987; Calvet et al. 1991; Chiang & Goldreich 1997). Fig¬ 
ure 9 illustrates some examples of the relationship be¬ 
tween the vertical distributions of p and T. 

Although it is a simplification, Eq. (5) faithfully high¬ 
lights this fundamental point about disk structures: the 
densities and temperatures are physically coupled. The 
intrinsic complexity of this coupling is that these physical 
conditions are linked to the radiative transfer of energy 
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Fig. 9.— Vertical cuts at a fixed radius (r = 10 AU) of the density 
(top) and temperature ( bottom ) profiles for three illustrative disk 
structures. The simple vertically isothermal behavior described 
in Eq. (5) is shown as gray dash-dot curves. The nominal scale 
height (corresponding to the midplane temperature), H & 0.4 AU, 
is marked with a vertical gray line. A more realistic variant, where 
particles with a power-law size distribution are well-mixed with 
the gas and the vertical structure has been iterated using Eq. (4) 
and a continuum radiative transfer code is shown as dashed black 
curves. Finally, the blue and red curves represent particles with 
1 pm and 1 mm sizes (respectively) from a similar model where the 
well-mixed assumption is relaxed: the solids are distributed ac¬ 
cording to a size-dependent balance between turbulent mixing and 
hydrostatic support, following the prescription of Dubrulle et al. 
(1995) (see Sect. 6.1). All models use the same E and material 
properties, and employ the RADMC-3D code. 


through the disk, which itself both depends on the den¬ 
sities and sets the temperatures. This feedback between 
structure and radiative transfer is a generic feature. 

5.1.2. Observational Constraints 

The most common approach used to measure the verti¬ 
cal structure of a disk employs the infrared SED as a di¬ 
agnostic (cf., Fig. 3), and relies on modeling the connec¬ 
tion between the height and temperature of the effective 
surface layer (the infrared “photosphere”; e.g., Kenyon 
& Hartmann 1987; Chiang & Goldreich 1997; D’Alessio 
et al. 1998). Since more energy is absorbed for a disk 
with a larger vertical extent, the temperature of the sur¬ 
face layer - and therefore the luminosity - scales with its 
characteristic height. The rough A r mapping of the 
infrared SED (see Sec. 3.1) means that the spectral slope 
reflects the radial variation of the surface layer height: 
more flaring makes a redder SED. Some examples of this 
effect are shown in Figure 10. However, although the 
SED does depend on the vertical structure, it alone can¬ 
not be used to unambiguously quantify any associated 
physical parameters (e.g., H): again, unresolved tracers 
do not robustly measure spatial properties. 
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Fig. 10.— A demonstration of the impact that the characteristic 
height has on key observables, including scattered light (top) and 
the SED (bottom), using a simple structure model and the RADMC-3D 
code. The models have structures described by Eq. (5), with iden¬ 
tical parameters aside from the shape of the radial H profile. Each 
H(r) is normalized at 1 AU (to ~0.04AU), and has a power-law 
scaling H(r) oc r' ,/J with if) = 1.10, 1.15, and 1.20. The most flared 
model (if) = 1.20) exhibits brighter mid-infrared emission and has 
a more prominent scattered light disk. 


The most valuable complement is a resolved map of 
scattered light emission (Sect. 3.2). This emission gen¬ 
erally appears as a bipolar pair of conical nebulae (rep¬ 
resenting the scattering layers on each side of the disk) 
with a dark waist (the disk interior); some representa¬ 
tive examples are shown in Figure 11. The curvature, 
or opening angle, of these nebulae is set by the shape 
of the scattering surface and thereby the flaring geome¬ 
try of the disk structure (e.g., Bastien & Menard 1990; 
Lazareff et al. 1990; Whitney & Hartmann 1992; Wood 
et al. 1998; D’Alessio et al. 1999; Takami et al. 2014). 
All else being equal, more flaring produces more curva¬ 
ture, larger nebulae, and a higher overall luminosity and 
polarization fraction for the scattered light. The images 
in Figure 11 illustrate some of these effects. 

In practice, the inference of vertical structure parame¬ 
ters from such data relies on some assumptions (or exter¬ 
nal constraints) about the viewing geometry and optical 
properties of the scatterers. For intermediate inclination 
angles, these degeneracies are usually explored in a re¬ 
stricted forward-modeling process (Roddier et al. 1996; 
Close et al. 1998; Grady et al. 1999, 2000, 2007; Potter 
et al. 2000; Augereau et al. 2001; Krist et al. 2002, 2005; 
Schneider et al. 2003; Pinte et al. 2008; Kusakabe et al. 
2012). Such ambiguities are diminished for extreme in¬ 
clinations: the variation of the surface layer height with 
radius can be determined most directly from the neb¬ 
ular curvature in the edge-on case (e.g., Burrows et al. 
1996; Lucas & Roche 1997; Krist et al. 1998; Padgett 
et al. 1999; Stapelfeldt et al. 1998, 2003; Cotera et al. 
2001; Wolf et al. 2003; Perrin et al. 2006; Glauser et al. 
2008), or more indirectly from the morphology of the ra¬ 
dial brightness profile for face-on disks (e.g., Krist et al. 
2000; Trilling et al. 2001; Debes et al. 2013). 

Using Eq. (5) as a parametric structure, models of re- 
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Fig. 11. — Examples of scattered light morphologies: ( clock¬ 
wise, from top left) optical/near-infrared HST images of GM Aur 
(Schneider et al. 2003), IM Lup (Pinte et al. 2008), HV Tau C 
(Stapelfeldt et al. 2003), and HH 30 IRS (Burrows et al. 1996; Wat¬ 
son Sz Stapelfeldt 2004, 2007; image courtesy of K. R. Stapelfeldt). 
These images are reproduced with permission from the authors and 
the AAS (© AAS) and Astronomy & Astrophysics (© ESO). 


solved scattered light data indicate H/r ~ 0.05-0.25, 
with only a modest radial variation (oc r 0 - 1-0 - 3 ). Sim¬ 
ilar estimates are made from the infrared SED. As a 
point of reference, the hydrostatic solution for a vertically 
isothermal disk predicts that H/r oc r 1 / 4 . It is worth a 
reminder that these are not actually measurements of 
the pressure scale height H , but rather inferences based 
on constraints of the characteristic height of the surface 
layer where dust scatters or absorbs starlight (which are 
not necessarily co-located). Although these are expected 
to roughly track one another, there are some additional 
effects that impact their relationship (Sect. 6.1). 

One such complication occurs when the characteris¬ 
tic height does not rise monotonically with r. If mate¬ 
rial closer to the star is distributed to sufficiently large 
heights, it can attenuate the radiation field at larger 
radii. This shadowing is usually attributed to a “puffed- 
up” inner disk edge (or wall), where internal energy gen¬ 
erated by viscous dissipation (for disks with high accre¬ 
tion rates; e.g., D’Alessio et al. 1998) or direct irradiation 
of the inner disk edge (Natta et al. 2001; Dullemond et al. 
2001; D’Alessio et al. 2005; Calvet et al. 2005) heats the 
material on sub-AU scales and inflates the local scale 
height. However, it may also be produced if vertical 
mixing is more vigorous in the inner disk (e.g., Dulle¬ 
mond & Dominik 2004a) or in the presence of small-scale 
substructures. Regardless of its origin, such obscuration 
makes the outer disk cooler, lowers the effective surface 
layer, and therefore reduces the infrared SED (particu¬ 
larly at longer A; a vertically enhanced inner edge actu¬ 
ally produces more near-infrared emission) and scattered 
light luminosity (Dullemond et al. 2001; Dullemond & 
Dominik 2004a). These effects have been inferred from 
SED morphologies (e.g., Meeus et al. 2001; Natta et al. 


2001; Dominik et al. 2003; Acke et al. 2004). A more 
striking, direct confirmation is found in scattered light 
images that show the flared outer disk emerging out of a 
shadow at large r (e.g., Garufi et al. 2014). 

Synoptic infrared photometry surveys are suggesting 
that small-scale vertical substructures may be common. 
Variability at some level is ubiquitous (see Cody et al. 
2014), but the timescales observed in the mid-infrared 
are surprising: they are much shorter than expected for 
the disk regions being traced (~few AU; see Fig. 3; Sitko 
et al. 2008; Morales-Calderon et al. 2009; Flaherty et al. 
2012; Kospal et al. 2012; Rebull et al. 2014). When 
near-infrared monitoring is also available, a “see-sawing” 
wavelength dependence has been observed, where the 
blue and red ends of the spectrum change in the oppo¬ 
site sense (Muzerollc et al. 2009; Espaillat et al. 2011). 
In rare cases this variability can be mapped in multi¬ 
epoch resolved scattered light observations (Stapelfeldt 
et al. 1999; Watson & Stapelfeldt 2007; Wisniewski et al. 
2008). Models of the amplitudes, timescales, and wave¬ 
length dependence of the variability suggest an origin at 
sub-AU radii, where modulations of the vertical structure 
(due to quickly evolving non-axisymmetric features, tur¬ 
bulence, or warps) shadow the disk at larger radii (e.g., 
Flaherty & Muzerolle 2010; Flaherty et al. 2013). 


5.2. Radial Structure 
5.2.1. Physical Overview 

The radial variation of the density structure is encapsu¬ 
lated in the surface density profile E (= f p dz). As was 
already noted, E is the most informative, essential as¬ 
pect of disk structure. The shape of E controls not only 
the likelihood of planet formation, but also the initial 
orbits and subsequent migrations of any planets during 
their formation epoch (e.g., Kokubo & Ida 2002; Ray¬ 
mond et al. 2005; Miguel et al. 2011). Moreover, E offers 
an indirect snapshot of the mass and (angular) momen¬ 
tum flows that underpin the global structural evolution 
of disk material (e.g., Hartmann et al. 1998). 

The standard theoretical models for E(r) are again 
based on the gas physics (presuming that the solids will 
follow the gas distribution; but see Sect. 6.1), using a 
balance of viscous and gravitational torques (Lin & Pa- 
paloizou 1980; Ruden & Lin 1986; Ruden & Pollack 1991; 
Stepinski 1998). For a thin Keplerian disk, E (and its 
evolution) can be determined by solving the viscous dif¬ 
fusion equation (Lynden-Bell & Pringle 1974) 


81 1 
dt 


3d_ 
r dr 


(v^£) 


(7) 


where v is the viscosity. There is a vast literature as¬ 
sociated with accretion flows and angular momentum 
transport in disks that centers around solving Eq. (7) 
for a variety of viscosity distributions (see the reviews by 
Pringle 1981; Papaloizou & Lin 1995; Turner et al. 2014). 
The typical assumptions, at least in studies making con¬ 
nections to observational data (although see Hueso & 
Guillot 2005), are that mass loss occurs only through 
accretion onto the star and that the viscosity is a (time- 
independent) radial power-law, v oc r 7 . The latter is 
motivated by the Shakura & Sunyaev (1973) prescrip¬ 
tion for turbulent viscosities in accretion disks, where 
v oc Cg/S2. In this case there are analytic similarity so- 
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Fig. 12.— {top): The surface density profiles inferred from 
mm/radio interferometric observations of nearby disks (Andrews 
et al. 2009, 2010b; Isella et al. 2009, 2010a; Guilloteau et al. 2011), 
assuming a spatially homogeneous opacity spectrum, a 1% dust- 
to-gas ratio, and a similarity solution structure (cf., Eq. 8). These 
profiles are broken at spatial scales that are not directly resolved 
(i.e., where E has been extrapolated). Blue boxes mark the stan¬ 
dard Weidenschilling (1977a) construction of the MMSN around 
the giant planets; the dashed black curve shows the normaliza¬ 
tion and shape advocated by Hayashi (1981). (bottom): The same 
information, but from the perspective of the cumulative mass dis¬ 
tribution (the integrated mass inside some radius r). 

lutions to Eq. (7), with a E profile that behaves like a 
power-law with an exponential taper at large radii, 



The characteristic scaling radius, r c , and the normaliza¬ 
tion depend on time and the initial conditions (e.g., size, 
mass, and mass flow rate; Lynden-Bcll & Pringle 1974; 
Lin & Pringle 1987; Hartmann et al. 1998). 

Given the uncertainty in these accretion disk models, 
empirical approximations are often preferred. The most 
common is an old idea based on a simplistic reconstruc¬ 
tion of E in the solar nebula, where the current planet 
masses are augmented back to solar composition and 
then smeared into annuli (e.g., Young 1901; Edgeworth 
1949; Kuiper 1956; Cameron 1962; Alfvcn & Arrhenius 
1970; Kusaka et al. 1970). For such a minimum mass so¬ 
lar nebula (MMSN) disk, E is simply a truncated power- 
law: the canonical Weidenschilling (1977a) construction 
of the MMSN is shown in Figure 12. Although there 
are fundamental conceptual flaws with the MMSN 8 , it 

8 Primary among these is the assumption of a static nebula and 


remains a standard benchmark in the field. 

5.2.2. Observational Constraints 

The most direct constraints on surface density profiles 
come from resolved (interferometric) observations of the 
(sub-)millimeter continuum emitted by the disk solids 
(cf., Sect. 3.1). Such emission typically has low optical 
depths, meaning the surface brightness scales with the 
product Kj, B,j(T)Y,. The standard approach to inter¬ 
pret the data is to assume that k„ is spatially homoge¬ 
neous and then fit or solve for T (as a radial power-law or 
with a radiative transfer calculation, respectively), and 
fit for E (as in Eq. 8 or a power-law), a characteristic or 
cutoff radius, and any other geometric parameters (e.g., 
inclination) with reference to the interferometric visibil¬ 
ities (e.g., Keene & Masson 1990; Lay et al. 1994, 1997; 
Wilner et al. 1996, 2000, 2003; Looney et al. 2000). For 
marginally resolved data, there are strong degeneracies 
between the radius parameter, inclination, and the radial 
gradients in T and E (e.g., Mundy et al. 1996). These 
can be mitigated with external constraints (the SED or 
scattered light data) on T and/or geometry (e.g., Akeson 
et al. 2002; Kitamura et al. 2002; Pietu et al. 2006; Hami- 
douche et al. 2006; Andrews & Williams 2007b; Pinte 
et al. 2008), but they persist at problematic levels unless 
the brightness distribution is resolved well. 

For the nearest disk populations, a resolution of <0'/5 
(60-70 AU) is usually sufficient to alleviate these degen¬ 
eracies and place rudimentary constraints on E. Surface 
density profiles have been inferred with such data for a 
few dozen disks, primarily in the Taurus (Isella et al. 
2009, 2010a; Pietu et al. 2006, 2014; Guilloteau et al. 
2011) and Ophiuchus (Andrews et al. 2009, 2010 b) re¬ 
gions (Fig. 12). With access to radii >15-30AU, these 
studies find a wide range of E gradients (7 ~ 0 - 1 ) and 
characteristic sizes ( r c ss 5-200 AU) 9 ; some of that diver¬ 
sity reflects the different modeling approaches. Figure 12 
also shows an alternative view in the form of cumulative 
mass distributions, which are perhaps more relevant for 
planet formation models. Many disks have roughly suffi¬ 
cient masses at 5-30 AU (based on model extrapolations) 
to form planetary systems like our own. However, the 
samples observed so far are biased (by necessity) toward 
the brightest, and thereby most massive, disks. 

Andrews et al. (2010b) noticed an interesting trend, 
now confirmed by Pietu et al. (2014), that fainter disks 
are typically more compact; see Figure 13. This crude 
size-luminosity (mass) relationship could be a manifes¬ 
tation of initial conditions (diversity in core angular mo¬ 
menta) or viscous timescales (Andrews et al. 2010 b); or 
perhaps it reflects variety in the evolutionary states of 
the solids (Pietu et al. 2014; see Sect. 6.1 for further dis¬ 
cussion). Following up on this trend is an active area of 

planetary system, meaning the MMSN density profile is a time- 
integrated structure (a convolution of two unknown effects: the 
density evolution and the efficiency by which material is accreted 
into the planets) rather than an instantaneous snapshot like would 
be observed in real disks or properly applied to a planet formation 
model (see also Cameron 1988; Davis 2005; Desch 2007). 

9 There is an often unappreciated (although obvious) difference 
between cutoff radii in power-law models of E and the character¬ 
istic radii ( r c ) in the similarity solutions (cf., Eq. 8): the latter 
encircle some significant fraction of the mass, but there is a non- 
negligible extension of low-density material at larger radii. Direct 
comparisons between different kinds of models require some care. 
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Fig. 13.— The 1.3 mm flux (scaled to a common distance), which 
scales with M as a function of the characteristic radius r c for all 
disks observed with a resolution better than 0775 and modeled with 
the similarity solution prescription for S (Eq. 8). Measurements are 
from Andrews et al. (2009, 2010b), Isella et al. (2010a), Guilloteau 
et al. (2011), and Pietu et al. (2014), and exclude binaries (see 
Sect. 4.2) and transition disks (see Sect. 6.2). 

research, but larger samples and ancillary observations 
will be required to test hypotheses for its origins. 

It is important to differentiate between what is being 
assumed and what is actually measured when interpret¬ 
ing resolved continuum data as described above. The sur¬ 
face brightness scales roughly with B v ( 1 — e~ Tu ) (Eq. 2), 
or ~ B V {T) t„ in the optically thin limit. A constraint on 
E requires a good (and internally self-consistent) model 
of the spatial variations of T (i.e., that the vertical struc¬ 
ture and irradiation source properties are well deter¬ 
mined) and opacity (recall t„ = E). For the latter, the 

typical assumption is that k v is spatially invariant, which 
is not especially well-justified (see Sect. 5.3). In fact, ra¬ 
dial and vertical variations in k v are generally expected 
(Sect. 6.1), and therefore will impact any estimates of 
E and T. While there are other caveats in this line of 
work, it is these basic systematic uncertainties about the 
material properties of the solids that ultimately limit the 
ability to robustly constrain disk E profiles. 

5 . 3 . Material Properties 

It should be obvious at this point that the links be¬ 
tween observations and the physical properties of disks 
are fundamentally tied to the detailed interactions of 
their constituent solid particles with radiation. Those 
interactions are affected by the material properties of 
the solids, particularly their compositions, sizes, and 
morphologies. This information is parametrically en¬ 
coded in the “optical” properties, absorption and scat¬ 
tering opacities as well as scattering and polarization 
phase functions, that control the redistribution of radia¬ 
tive energy in both frequency and direction. These can 
be determined for any collection of particle properties 
and their corresponding optical constants (wavelength- 
dependent refractive indices) using numerical methods 
of varying sophistication, from simplified schemes like 
Mie theory (e.g., Bohren & Huffman 1983) to the brute- 
force simulated solutions of Maxwell’s equations in the 
discrete dipole approximation (Purcell & Pennypacker 
1973; Draine & Flatau 1994). The optical constants for 
many relevant materials have been calculated theoreti¬ 


cally (e.g., Draine & Lee 1984; Laor & Draine 1993) or 
measured in the laboratory (e.g., Jaeger et al. 1994,1998; 
Dorschner et al. 1995; Henning et al. 1995). 

5 . 3 . 1 . Particle Composition 

Reasonable models for the composition of disk solids 
(e.g., Pollack et al. 1994) are based on assumed chemi¬ 
cal pathways and abundance measurements in primordial 
solar system bodies (Lodders 2003): the primary contrib¬ 
utors include silicates, carbonaceous materials, metallic 
compounds, and water ice. The emission bands at ~10 
and 18 /im in disk spectra verify an abundant popula¬ 
tion of small silicate particles with a range of miner¬ 
alogies (Forrest et al. 2004; Kessler-Silacci et al. 2005, 
2006; Watson et al. 2009; Sargent et al. 2006, 2009b,a; 
Olofsson et al. 2009; Oliveira et al. 2010, 2011). Addi¬ 
tional features throughout the infrared provide evidence 
for hydrocarbons (PAHs; e.g., Habart et al. 2004; Geers 
et al. 2006, 2007; Keller et al. 2008), water ice (Chiang 
et al. 2001; Creech-Eakman et al. 2002; McClure et al. 
2015), and other relevant minerals (e.g., Sturm et al. 
2013). Resolved maps of scattered light colors suggest 
organic grain coatings (e.g., tholins; Debes et al. 2008, 
2013; Rodigas et al. 2015). While it is not trivial to quan¬ 
tify the relative abundances of these materials, particu¬ 
larly in the midplane where spectral features are hidden 
by high optical depths, the combined constraints from 
disk spectra, analyses of meteorites, comets, and aster¬ 
oids, and chemical intuition provide a reasonably coher¬ 
ent picture of the bulk composition of disk solids. 

5 . 3 . 2 . Particle Size Distribution 

The particle size (radius), a, dictates the cross-section 
for interaction with radiation of a given wavelength, and 
thereby has a significant impact on k v . The opacity spec¬ 
trum is approximately constant for a A (geometric op¬ 
tics limit), falls off like A -2 when a <C A (Rayleigh limit), 
and is enhanced by resonant interactions at intermediate 
wavelengths. The implication is that emission observed 
at a given wavelength most efficiently probes particles 
of a similar size, a ~ A. In a realistic disk environment 
there will be a distribution of particle sizes, often approx¬ 
imated as a power-law, dN/da oc a~ q . Typically, most of 
the mass is in the largest particles (q < 4); as a reference 
point, simple collisional models (Dohnanyi 1969) and fits 
to the interstellar extinction curve (Mathis et al. 1977; 
Cardclli et al. 1989; Mathis 1990) suggest q ~ 3.5. 

Figure 14 illustrates how the opacity spectrum changes 
for different particle size distributions. As the maximum 
size increases (or the distribution becomes more top- 
heavy), there are three general changes to the opacities. 
First, the infrared opacities decrease. However, if the 
disk is still optically thick in the infrared, it is not easy 
to demonstrate this observationally. Second, the silicate 
spectral features become muted (e.g., Min et al. 2004). 
This has been identified as a signature of grain growth 
in the disk surface layers (e.g., Bouwman et al. 2001; 
Meeus et al. 2001; Kessler-Silacci et al. 2006; Olofsson 
et al. 2009; Lommen et al. 2010), although there is room 
for structural (e.g., Dominik & Dullemond 2008) or com¬ 
positional (Honda et al. 2003; Meeus et al. 2003) ambi¬ 
guity in quantifying the size distribution. And third, the 
approximately power-law spectral behavior at mm/radio 
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Fig. 14.— {top): Examples of absorption opacity spectra for sim¬ 
ple models of dust populations relevant for disks, based on the ho¬ 
mogenized analysis efforts of the DIANA project (Woitke et al. 2015). 
The composition is a mixture of amorphous silicates, carbon, and 
vacuum, with various assumptions about the size distribution and 
porosity (in this case, “compact” and “fluffy” mean filling factors 
of 1 and 0.5, respectively). Optical constants were collated from 
Dorschner et al. (1995) and Zubko et al. (1996), and combined us¬ 
ing the Bruggeman (1935) mixing rule. The ‘distribution of hollow 
spheres’ approach advocated by Min et al. (2005) was then used to 
compute opacities, (bottom): The 1-3 mm opacity spectral index, 
/3, as a function of the maximum particle size. The typical spectral 
slopes measured for disks cluster around /3 < 1, indicating that 
most of the particle mass has grown to at least mm sizes. 

wavelengths flattens; i.e., if K mm oc then /3 decreases 
with larger a max or smaller q (see Fig. 14; Miyake & 
Nakagawa 1993; Ossenkopf & Henning 1994; Henning & 
Stognienko 1996; D’Alessio et al. 2001; Draine 2006). 

This latter point is especially interesting for disks be¬ 
cause they become optically thin at these long wave¬ 
lengths. The spectral dependence of the mm/radio con¬ 
tinuum therefore scales with B v k„; in the Rayleigh- 
Jeans limit the observed spectrum would scale like v 2+ P , 
and so the radio “colors” serve as a key diagnostic of the 
particle size distribution. 10 For the small (sub-micron) 
grains in the interstellar medium, the mm/radio spec¬ 
trum is steep; /3 ~ 1.8±0.2 (Hildebrand 1983; Goldsmith 
et al. 1997; Finkbeiner et al. 1999; Li & Draine 2001). 
The integrated spectra of ~Myr-old disks have compar¬ 
atively much “redder” colors, indicative of (3 < 1 in most 
cases (Beckwith & Sargent 1991; Mannings & Emerson 
1994; Andrews & Williams 2005; Wilncr et al. 2005; Lorn- 
men et al. 2007, 2009; Ricci et al. 2010b,a; Ubach et al. 
2012 ) and therefore substantial growth (a max ~nun or 
cm; e.g., Draine 2006). There are a few important points 
to keep in mind regarding the interpretation of spectral 
behavior in this context. For shallow spectra (low /3), the 

10 In practice, the Rayleigh-Jeans limit is often not strictly valid 
for disks. However, even simple assumptions for T can be used to 
appropriately account for the curvature in the Planck term. 


relationship between the opacity spectral index and par¬ 
ticle size saturates (see Fig. 14); when solids reach sizes 
much larger than the observing wavelengths, they emit 
much less and have negligible impact on the observed 
spectrum. So, associating /3 with a particle size distri¬ 
bution extending beyond cm-sized solids is generally not 
possible. More practically, it is crucial to spatially resolve 
brighter disks, so any contamination from optically thick 
emission (with a spectrum that scales like v 2 ) can be 
identified (Testi et al. 2003; Natta et al. 2004; Rodmann 
et al. 2006; Ricci et al. 2012). Likewise, when increasing 
the spectral range to the radio (A ~ 1 cm) in an effort 
to better leverage constraints on /?, the contamination 
of non-disk emission (free-free or synchrotron radiation) 
also needs to be considered (e.g., Rodmann et al. 2006). 

5 . 3 . 3 . Particle Morphology 

The optical properties of solids are also sensitive to 
their structure or morphology, typically discussed in 
terms of porosity (or the volume filling factor; e.g., 
Kozasa et al. 1992; Henning & Stognienko 1996; Min 
et al. 2003; Kimura et al. 2003; Orrnel et al. 2011). Nu¬ 
merical simulations predict that small individual grains 
(monomers) will stick together in low-speed collisions 
and grow into fluffy, porous aggregates (e.g., Ossenkopf 
1993; Dominik & Tielens 1997; Ormel et al. 2007; Wada 
et al. 2009; Okuzumi et al. 2012); laboratory experiments 
directly confirm this behavior (see Blum & Wurm 2008). 
The optical properties of porous aggregates depend on 
the product of their bulk size and filling factor (e.g., 
Kataoka et al. 2014). The absorption opacities of large 
aggregates with high porosity are similar to those for 
smaller, compact particles, with the subtle difference be¬ 
ing that the resonant enhancement of for a ~ A is 
muted (cf., Fig. 14). The corresponding scattering opac¬ 
ities may be relatively enhanced at long wavelengths: 
Kataoka et al. (2014) suggest mm-wave polarimetry as 
a potential probe (although current limits are already 
quite strong; e.g., Hughes et al. 2009a, 2013). More infor¬ 
mation on the morphologies of disk solids are available, 
at least in the surface layers, through constraints on the 
scattering phase function (e.g., Mulders et al. 2013) and 
polarized intensities (e.g., Min et al. 2012). 

5 . 3 . 4 . Caveats on Mass and Density Estimates 

On consideration of all these complex microphysical 
factors that set the opacities, which in turn are respon¬ 
sible for translating physical conditions into the astro¬ 
nomical tracers that are observed, the effort to quantify 
the role of disks in planet formation may seem incredibly 
daunting. And indeed, there is much work to be done in 
teasing out the complex inter-dependences between the 
material and structural properties of disks. 11 So, it is 
imperative to consider inferences about disk properties 
in the context of the assumed (i.e., mo del-dependent) 
properties of their constituent solids. 

Given the fundamental significance of disk masses in 
the planet formation process, there is special interest in 
understanding the uncertainty (or diversity) in the opac¬ 
ities at millimeter wavelengths. With the few exceptions 

11 The situation is in fact more complicated than presented, since 
the compositions, sizes, and morphologies of the solids all vary as 
a function of the physical conditions (e.g., see Sect. 6.1). 
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of disks with especially simple structures (e.g., Andrews 
et al. 2014), this issue is best assessed demographically. 
A few different arguments have been made that point 
toward the “standard” millimeter opacities (~2 cm 2 g -1 
at 1.3 mm; Beckwith et al. 1990) being over-estimated, 
and therefore that the disk masses may be higher than 
are typically inferred. One points out the order of mag¬ 
nitude mismatch between the masses inferred from the 
millimeter emission and a crude proxy using the prod¬ 
uct of accretion rates and ages (Hartmann et al. 1998, 
2006). Another derives from the claim that the massive 
disk fraction seems to be too low to produce the observed 
giant planet frequency (e.g., Greaves & Rice 2010; Na- 
jita & Kenyon 2014). 12 The direction of this proposed 
uncertainty makes sense, since the growth of solids will 
lock up a fraction of the mass in large bodies that do 
not emit much at millimeter wavelengths (cf., Fig. 14). 
But its scale is perhaps less clear. It is unlikely that 
the masses could be substantially higher than estimated 
without a large population of gravitationally unstable 
disks, for which no evidence (e.g., Narayanan et al. 2006; 
Jang-Condell & Boss 2007; Dipierro et al. 2014) has yet 
been identified (see especially the symmetry for the very 
bright HL Tau disk; ALMA Partnership et al. 2015). So 
considering these (admittedly approximate) constraints, 
it seems reasonable to assume that the systematic uncer¬ 
tainty on the mean millimeter opacity in disks (relative to 
the typically adopted value) is modest, about (±)ldex. 

5 . 4 . Synopsis 

• The flared, vertical distribution of dust grains reg¬ 
ulates the thermal structure of a disk. The broad¬ 
band infrared SED and spatially resolved measure¬ 
ments of optical/infrared scattered light can be 
used to constrain the vertical structure. In the con¬ 
text of simple models, such data suggest that the 
typical vertical aspect ratio (~U/r) is 0.1. 

• The radial surface density profile, E, is a crucial 
input for planet formation models. Interferomet¬ 
ric measurements that resolve the optically thin 
millimeter continuum emission can be used to in¬ 
terpret the radial distribution of solid mass; some 
preliminary results are consistent with models of 
the (outer) primordial solar nebula. 

• The compositions, sizes, and shapes of the solid 
particles determine the opacities, which regulate 
the radiative transfer of energy into and out of 
the disk structure. There are strong and com¬ 
plex inter-dependencies between these material and 
structural properties. Measurements of the shape 
of the mm/radio continuum spectrum are sensi¬ 
tive to the particle size distribution; the relatively 
“red” mm/radio colors that are observed indicate 
the substantial growth of disk solids. 

• The physical conditions inferred from the interpre¬ 
tation of disk observations need to be considered 
only in the context of the model assumptions (e.g., 

12 This discrepancy is not too large, since roughly a quarter of 
~Myr-old disks around Sun-like hosts have masses in excess of the 
MMSN (Andrews et al. 2013; Fig. 6). Nevertheless, it demonstrates 
that the masses are likely not systematically over-estimated. 


the functional form of the density distribution, the 
material properties of the constituent solids). 

Additional Reading', the review of disk modeling by 
Dullemond et al. (2007); a basic review of resolved struc¬ 
ture constraints by Wilner & Lay (2000); the more gen¬ 
eral review by Williams & Cieza (2011). 

6. EVOLUTION 

The general framework for interpreting protoplanetary 
disk observations, as outlined in the previous sections, is 
still rather crude and over-simplified. With some recent 
enhancements in data quality, it has become clear that 
a more nuanced approach is necessary. There is now un¬ 
equivocal evidence that evolutionary effects play a large 
(and likely dominant) role in determining the observable 
characteristics of these disks. The goal of this section is 
to highlight how two key evolutionary behaviors - the 
growth and migration of solids, and the relatively rapid 
clearing of the inner disk impact the standard obser¬ 
vational benchmarks of disk structure, and thereby to 
point out opportunities for new insights on the physical 
mechanisms that drive such evolution. 

6 .1. The Growth and Migration of Solids 

Most observational studies of disk structures implicitly 
assume that the collisional evolution of the constituent 
solids can be ignored. That is a bad assumption. The 
growth and migration of these solids involve extraor¬ 
dinarily complex processes, but exploring their conse¬ 
quences - both physically and observationally - is a cru¬ 
cial step toward developing a more comprehensive un¬ 
derstanding of planet formation. To illustrate the effects 
of this evolution, in this section we will consider what 
happens to an ‘initial’ population of small dust grains 
suspended in a (hydrostatically supported) gas reservoir 
with a homogeneous dust-to-gas ratio, £(r, z) = Co- 

6 . 1 . 1 . Settling 

Relative motions imparted stochastically (Brownian 
motion, turbulent mixing) will make these grains col¬ 
lide; at low speeds, they will stick together and produce 
larger aggregates (e.g., Blum & Wurrn 2008; Guttler et al. 
2010). Smaller particles are coupled to the gas, buoyed 
up into the disk atmosphere by hydrostatic support. But 
larger particles feel less of that lift and will settle toward 
the midplane (Adachi et al. 1976; Weidenschilling 1977a; 
Nakagawa et al. 1981). As they sink into denser regions, 
the collision rate increases and growth can accelerate. 
This growth and (vertical) migration feedback process 
is not perfectly efficient (Dullemond & Dominik 2005): 
alternative collision outcomes like electrostatic repulsion 
(Okuzumi 2009) and bouncing (Zsom et al. 2010) can 
slow growth, while fragmentation (Brauer et al. 2008; 
Birnsticl et al. 2009) and erosion (e.g., Seizinger et al. 
2013; Krijt et al. 2015) can actually reverse it. 

The interactions between the gas pressure, gravity, and 
(turbulent) mixing, coupled with the efficiency of parti¬ 
cle growth (relative to destruction), gives rise to a size- 
sorted, sedimentary layering of the disk solids in which 
( decreases with z (e.g., Volk et al. 1980; Cuzzi et al. 
1993; Dubrulle et al. 1995; Schrapler & Henning 2004). 
This size segregation increases the vertical temperature 
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gradient and concentrates the absorption/scattering op¬ 
tical depth profile toward the midplane, reducing the ir¬ 
radiated surface area and height of the surface layer and 
thereby impacting some key observables (e.g., Dullemond 
& Dominik 2004b; see Sect. 5.1). 

A settled disk has a lower characteristic surface layer, 
absorbs less starlight, and therefore is colder. Conse¬ 
quently, it emits less thermal and scattered continuum 
radiation to the observer (Miyake & Nakagawa 1995; 
Dullemond & Dominik 2004b; D’Alessio et al. 2006; Mul¬ 
ders & Dominik 2012). Distinguishing settling from a 
disk with an intrinsically low gas pressure scale height, 
H , is not trivial. Some indirect constraints can be made 
from comparisons of mid-infrared colors and the strength 
of the 10 fim silicate feature (e.g., Kessler-Silacci et al. 
2006; Furlan et al. 2009), or the inconsistent inferences 
of surface height layers from different tracers (e.g., scat¬ 
tered light images and infrared SEDs; Stapelfeldt et al. 
2003; Wolf et al. 2003). A more direct signature of this 
stratification is accessible by comparing scattered light 
images at multiple wavelengths. Since dust grains pref¬ 
erentially scatter at wavelengths comparable to the par¬ 
ticle size, settled disks show a scattering layer height that 
decreases with A (Duchene et al. 2003, 2004, 2010; Pinte 
et al. 2007; Grafe et al. 2013). Edge-on disks may be 
particularly useful in future work on this subject, when 
ALMA can spatially resolve both the temperature struc¬ 
ture of the gas (e.g., Dartois et al. 2003; Rosenfeld et al. 
2013a; de Gregorio-Monsalvo et al. 2013) and the contin¬ 
uum that traces ~mm particles as a function of height 
above the midplane (e.g., Boehler et al. 2013). As a ref¬ 
erence, radiative transfer models of typical disk SEDs 
indicate that settling can reduce the dust-to-gas ratio in 
the disk atmosphere by a factor 1/e ~ 10-1000, where 
C = e (o in the parlance of D’Alessio et al. (2006). 

6.1.2. Radial Drift 

Similar evolution also proceeds in the radial dimension. 
The radial motions of disk solids are fundamental to that 
evolution. Particles that are well coupled to the gas will 
follow its viscous flow trajectories, moving inward toward 
the host star at small radii (accretion), and away from 
it in the outer disk (diffusion) (e.g., see Takeuchi & Lin 
2002). These radial velocities induced by viscous drag 
tend to be slow (~few to tens of cm s -1 ), and are most 
relevant for small (~^m-sized) particles. An aerody¬ 
namic process termed “radial drift” can be substantially 
more influential. In general, the gas motions in a disk 
are dominated by rotation around the host star. How¬ 
ever, the presence of a radial pressure gradient imparts a 
small deviation from Keplerian velocities. For example, 
in the standard picture of a smooth disk with temper¬ 
atures and densities that decrease with r, the pressure 
gradient is negative and the gas disk rotates at a slightly 
sub-Keplerian rate. In the limit of weak coupling to the 
gas, solids orbit at Keplerian speeds. In the intermedi¬ 
ate drag regime, particles with a certain range of sizes 
feel a strong headwind due to their velocity differential 
with the gas, which saps their orbital energy and sends 
them spiraling toward the pressure maximum (Whipple 
1972; Adachi et al. 1976; Weidenschilling 1977a; Naka¬ 
gawa et al. 1981, 1986). This radial drift moves solids so 
efficiently that they can be locally depleted faster than 
the coagulation timescale: in a MMSN disk, such a ‘drift 



Fig. 15.— Demonstration of the evolution of solids embedded 
in a gas-rich disk. In each panel, grays refer to the gas, blues to 
~/mi-sized grains, and reds to ~mm-sized particles. Three time 
steps are shown - 0.5, 1, and 2Myr (from light to dark) - start¬ 
ing from an initial homogeneous distribution of sub-/mi monomers 
with £ = 0.01 in a standard viscous accretion disk (see Birnstiel & 
Andrews 2014 for details), (top) Radial velocities: large negative 
(inward) motions (up to ~0.1 km s' -1 levels) due to radial drift are 
noted for mm-sized solids at tens of AU. ( middle ) Surface densities 
(E): the mm solids are radially concentrated and have a relatively 
sharp “edge” well interior to the gas or small grain distributions. 
(bottom) Dust-to-gas mass ratio (£): a steep drop in £ m m is pro¬ 
duced by radial drift. But note the scale in these low values are 
emblematic of a major issue with the radial drift timescales (they 
are much too short) in current evolutionary models. 


barrier’ to further growth occurs at ^1AU for m-sized 
bodies, and at ~100AU for mm-sized particles. 

The radial migration outlined above significantly mod¬ 
ifies the distribution of disk solids in both size and 
space. Figure 15 illustrates these effects for a representa¬ 
tive evolutionary model (cf., Birnstiel & Andrews 2014). 
The higher densities and relative velocities present at 
smaller disk radii naturally enhance the local collision 
(and thereby growth) rates (e.g., Dullemond & Dominik 
2005). But the added influx of ~mm/cm-sized particles 
from the outer disk will amplify this radial size-sorting. 
The net result is that larger solids will be preferentially 
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Fig. 16.— Illustration of the observational effects of particle 
growth and migration in the TW Hya disk, (top, left to right) 
The CO J= 3—2 line (color scale is the intensity-weighted velocity 
field), 870 continuum, and 9mm continuum emission observed 
by the SMA (Andrews et al. 2012) and VLA (Menu et al. 2014) 
interferometers. The images have the same scale (4", ~430AU, 
on a side) and resolution (l", ~50AU). Note the progression of 
radial concentration in the various tracers that is predicted by the¬ 
oretical models (cf., Fig. 15). (bottom) The same progression from 
extended (CO) to compact (9mm continuum), now viewed more 
directly as the azimuthally averaged visibility profiles (see Sect. 3.1 
and Fig. 4). The low-level oscillations in the 870 visibilities are 
a distinct signature of the “sharp” edge in the radial distribution of 
~mm-sized particles imposed by radial drift (see Fig. 15, bottom , 
and Birnstiel & Andrews 2014). 

concentrated close to the gas pressure maximum (typ¬ 
ically the inner disk edge; e.g., Takeuchi & Lin 2002; 
Brauer et al. 2008; Birnstiel et al. 2010, 2012). Since most 
of the solid mass is in larger particles that are system¬ 
atically depleted from the outer disk by radial drift, the 
dust-to-gas ratio C decreases with r (Takeuchi et al. 2005; 
Birnstiel & Andrews 2014). For the sizes where drift ve¬ 
locities are maximal in the outer disk (~mm/cm), the 
ratio (Cmm) drops precipitously at ~10-100AU (Youdin 
& Shu 2002; Weidenschilling 2003; Takeuchi & Lin 2005; 
Jacquet et al. 2012; Birnstiel & Andrews 2014). 

6.1.3. Signatures of Radial Drift 

These physical effects have pronounced observational 
hallmarks. The radial sorting of particle sizes induced by 
growth and drift is manifested in the disk opacities, 
particularly at mm/cm wavelengths (cf., Sect. 5.3). Since 
larger particles have a lower opacity spectral index, /3 
(a “redder” n v ), disks that have experienced substantial 
evolution in their solids will have an increasing f3(r): such 
an opacity profile is observed as a radial variation of the 
radio colors (i.e., a relatively “blue”, or steep, spectrum 
at larger r; Isella et al. 2010a). In practice, this color 
gradient is observable as an anti-correlation between the 
wavelength and size of the continuum emission, such that 
the intensity distribution at longer wavelengths is more 
compact (Banzatti et al. 2011; Guilloteau et al. 2011; 
Perez et al. 2012, 2015; Trotta et al. 2013; Menu et al. 
2014). Figure 16 shows an example of this effect. 

A confirmation of the related prediction that drift pro¬ 


duces a decreasing ((r) (cf., Fig. 15) is more difficult, 
since accessing and interpreting gas density tracers re¬ 
mains a challenge. However, some indirect, qualitative 
evidence supporting a steep drop-off in Cmm (r) is avail¬ 
able through a comparison of the mm/cm-wavelength 
continuum and molecular line emission sizes: the for¬ 
mer should have a more radially compact distribution 
than the latter. Such line/continuum size discrepancies 
have been noted for some time (e.g., Pietu et al. 2005; 
Isella et al. 2007), but were (reasonably) attributed to op¬ 
tical depth effects: since the mm/cm continuum is opti¬ 
cally thin, insufficient sensitivity can lead to a size under¬ 
estimate relative to the optically thick lines that typically 
trace the gas (Hughes et al. 2008). As the data qual¬ 
ity has improved, it has nevertheless become clear that 
these apparent size discrepancies are not observational 
artifacts, but are instead intrinsic to the disk structures. 
For a still-limited sample, this signature of radial drift 
is manifested in a mm continuum distribution ~2-5x 
smaller than for the CO line (Panic et al. 2009; Andrews 
et al. 2012; Rosenfeld et al. 2013b; de Gregorio-Monsalvo 
et al. 2013; Walsh et al. 2014; see Fig. 16). 

Some complementary information is encoded in the 
morphology of the mm continuum at larger radii. The 
steep drop in ( mm (r) predicted theoretically to be a 
consequence of radial drift (e.g., Birnstiel & Andrews 
2014; see Fig. 15) should produce a relatively sharp 
outer “edge” in the continuum distribution. That edge is 
clearly identifiable in the oscillatory structure it imprints 
on the interferometric visibilities (see Fig. 16; Andrews 
et al. 2012; de Gregorio-Monsalvo et al. 2013). 

6.1.4. Open Issues 

The dramatic modifications these evolutionary factors 
impart on the traditional observational tracers of struc¬ 
ture (cf., Sect. 5), should make it clear that constraining 
disk densities is a substantially more complicated task 
than has often been assumed. At some level, the spatial 
variations imposed on k v and £ by these evolutionary 
mechanisms render many of the standard assumptions 
used in the interpretation of disk structures invalid. To 
put it bluntly, there is not currently any obvious, coher¬ 
ent parametric framework available for characterizing the 
density distributions of disk solids. Right now, the focus 
is instead on collecting additional evidence that supports 
the evolution of solids, and building up a foundation for 
its proper interpretation. The expectation is that re¬ 
solved multiwavelength continuum and spectral line data 
can be combined to better characterize the growth and 
migration of solids, and thereby to facilitate integrated 
studies of structure and evolution. 

This subject is a relatively new, and incredibly dy¬ 
namic, emphasis for the held. But while the influx of 
high-quality resolved data is qualitatively in sync with 
theoretical predictions for the growth and migration of 
disk solids, it is important to recognize that there are 
also some serious conflicts. Most prominent among these 
is the issue of migration timescales. Theory predicts that 
radial drift is much too efficient to explain the spatial dis¬ 
tributions and luminosities of the mm-wavelength contin¬ 
uum that are routinely observed for ~Myr-old disks (e.g., 
Takeuchi & Lin 2005; Brauer et al. 2007). The hypoth¬ 
esized remedy to relieve this tension is that there must 
be small-scale substructure in the gas disk (e.g., Whipple 
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1972; Pinilla et al. 2012b). Localized pressure maxima 
introduced by turbulent fluctuations (Klahr & Henning 
1997; Dzyurkevich et al. 2013; Flock et al. 2015), mag¬ 
netohydrodynamic (MHD) structures (i.e., zonal flows; 
Johansen et al. 2009; Bai & Stone 2014; Suzuki & Inut- 
suka 2014; Bai 2015), condensation fronts (e.g., Kretke & 
Lin 2007), or dynamical perturbations from (planetary) 
companions (Pinilla et al. 2012a; Birnstiel et al. 2013) 
can slow or even “trap” drifting particles, potentially 
reconciling the theoretical timescales and observational 
constraints. There is some preliminary evidence in sup¬ 
port of these kinds of substructure for a special subset of 
disks (see Sect. 6.2), and there will soon be data capa¬ 
ble of quantifying its nature and prevalence in the more 
general disk population (see Sect. 7). 

6 .2. The “Transition” Phase 
6.2.1. Physical Overview 

The evolution of angular momentum, through turbu¬ 
lent transport (Shakura & Sunyaev 1973; Lynden-Bell 
& Pringle 1974; Hartmann et al. 1998) or dissipation 
via magnetic winds (Blandford & Payne 1982; Konigl 
1989; Bai & Stone 2013), largely dictates the long-term 
behavior of disk structures. If acting alone, such pro¬ 
cesses would steadily deplete the densities until the disk 
slowly (tens of Myr) faded away. However, there are ad¬ 
ditional mechanisms that should expedite the dissipation 
or transformation of the disk material. 

In terms of dispersal, the primary pathway of interest 
is through a photoevaporative wind. Eneregtic radiation 
(far-ultraviolet to X-ray) generated by the host star can 
impart the gas in the disk surface with enough energy to 
escape the system (Clarke et al. 2001; Alexander et al. 
2006a; Ercolano et al. 2008; Gorti & Hollenbach 2009; 
Owen et al. 2010). When the mass-loss rate in the re¬ 
sulting wind exceeds the (inward) viscous flow rate, a 
gap opens near the wind launch radius (^few AU; e.g., 
Liffman 2003; Font et al. 2004; Ercolano et al. 2009). 
The inner disk is cut off from re-supply, and will quickly 
(<0.1 Myr) accrete onto the star (~0.1 Myr; e.g., Alexan¬ 
der et al. 2006b; Gorti et al. 2009; Owen et al. 2011). This 
depletion makes the inner disk optically thin in the radial 
direction, permitting a more efficient irradiation at the 
inner edge of the remaining disk. This in turn enhances 
the wind mass-loss rate and thereby quickly (<0.01 Myr) 
expands the wind-driven “cavity” at the disk center to 
r > 10 AU (Clarke et al. 2001; Alexander et al. 2006b; 
Owen et al. 2012); it may even trigger an instability that 
destroys the disk altogether (Owen et al. 2013). 

In the alternative (but not mutually exclusive; e.g., 
Matsuyama et al. 2003; Rosotti et al. 2013) context of 
metamorphosis, planet formation can also dramatically 
alter the disk structure on relatively short timescales. 
Once a planet is sufficiently massive (depending on 
the disk viscosity, but typically ^0.1-1 Mj up ), it exerts 
torques that repel the local disk material away from its 
orbital radius (e.g., Goldreich & Tremaine 1980; Lin & 
Papaloizou 1979, 1993; Bryden et al. 1999; Kley & Nel¬ 
son 2012). The net result is a narrow annulus of depleted 
densities, a gap. Depending on the local disk properties 
and planet mass, some disk material may continue to 
flow to the planet (or beyond it) in accretion streams 
(Artymowicz & Lubow 1996; Kley 1999; Lubow et al. 


1999). A more massive planet accretes a larger fraction 
of this flow (e.g., Lubow & D’Angelo 2006), diminishing 
the re-supply of the disk region interior to the gap, and 
thereby effectively creating a “cavity” at the disk cen¬ 
ter (e.g., Varniere et al. 2006; Crida & Morbidelli 2007; 
Dodson-Robinson & Salyk 2011; Zhu et al. 2011). 

Theoretical models of these processes make some im¬ 
portant and testable observational predictions. The most 
fundamental of these is that disk evolution should involve 
two timescales — a long stage where the processes related 
to momentum transport dominate, followed by a rapid 
“transition” phase where the disk structure is severely 
modified by interactions with winds and/or (giant) plan¬ 
ets — that can be identifiable in simple demographic 
data (e.g., see Alexander & Armitage 2009). The disks in 
this latter evolutionary phase have a distinct structural 
feature: a central region of depleted densities. Such a 
structure necessarily imposes a pressure maximum just 
outside this low-density cavity, which may concentrate 
drifting solids in a ring-shaped “trap” (cf., Sect. 6.1; e.g., 
Alexander & Armitage 2007; Pinilla et al. 2012a; Gorti 
et al. 2015). As with normal disks (cf., Sect. 6.1), these 
solids leave behind a radially extended gas disk with a 
low dust-to-gas ratio. The shape of the pressure gradi¬ 
ent around the maximum regulates the flow of material 
to smaller radii. If the gradient is not too steep, gas and 
smaller particles can pass through into the cavity; how¬ 
ever, larger solids are filtered out of this flow and remain 
trapped near the pressure maximum (e.g., Rice et al. 
2006; Paardekooper & Mcllema 2006; Zhu et al. 2012, 
2014). In a sense, the small amount of material in the 
cavity (or, rather, the flow rate of disk material through 
that pressure maximum), provides important clues to the 
physical mechanism(s) responsible for the depletion that 
is central to this transitional phase of evolution (e.g., 
Alexander et al. 2006b; Najita et al. 2007, 2015; Alexan¬ 
der & Armitage 2009; Owen & Clarke 2012). 

6.2.2. Observational Signatures 

The observational evidence for this rapid stage of disk 
evolution pre-dates (and motivates) much of the theoret¬ 
ical work. In any young cluster, there is a small sub¬ 
population (~1 10%) that lacks a near-infrared excess 
but still exhibits strong dust emission at longer wave¬ 
lengths (Strom et al. 1989; Skrutskie et al. 1990). These 
transition disks are usually identified by the “dip” in their 
infrared spectra, with blue colors shortward of ^5-10 pm 
and red colors at longer wavelengths (e.g., Brown et al. 
2007; Cieza et al. 2007, 2012; Merfn et al. 2010; Furlan 
et al. 2011). A representative example is shown in Fig¬ 
ure 17. Given the rough A H > r mapping between the 
spectrum and dust temperatures (Sect. 3.1), this “dip” 
signature suggests a substantial depletion of solids in the 
warm inner disk (Strom et al. 1989; Skrutskie et al. 1990; 
Marsh & Mahoney 1992, 1993; Calvet et al. 2002, 2005; 
Rice et al. 2003; D’Alessio et al. 2005). If all disks ex¬ 
perience such a transition, a simple duty cycle argument 
demands that it occurs quickly (the product of the tran¬ 
sition disk fraction and the mean sample age, so <0.1 
0.5 Myr; but see Sicilia-Aguilar et al. 2008; Muzerolle 
et al. 2010). Comparing that rough timescale with the 
decay time of infrared excesses in the general disk pop¬ 
ulation (Sect. 4.4; see Mamajek 2009), these simple in¬ 
frared color demographics clearly reveal the theoretically 
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Fig. 17.— The SED and representative resolved images of the 
transition disk around the young star SR21. The blue curve is a 
model of the stellar photosphere; the gray curve is a well-sampled 
spectrum from Spitzer. The SED shows a prominent “dip” in the 
infrared, produced by the substantial depletion of a cavity with 
radius ~20 AU. The image insets show the 1.6 //rn (polarized) scat¬ 
tered light (Follette et al. 2013), along with the 450 /im (Perez et al. 
2014) and 870 (im thermal continuum emission (Brown et al. 2009). 
Each image spans 175 AU on a side. The thermal images clearly 
show the low-density cavity, but there is scattered light tracing 
small grains located well inside the (sub-)mm ring. 

predicted two-timescale behavior for disk evolution. 

Moreover, the “dust ring” structures implied by these 
distinctive spectral morphologies have been directly con¬ 
firmed, and characterized in detail, using resolved ob¬ 
servations of the mm continuum (e.g., Pietu et al. 2005, 
2006; Hughes et al. 2007, 2009b; Brown et al. 2008, 2009, 
2012; Isella et al. 2010a,b; Andrews et al. 2009, 2010b, 
2011; Cieza et al. 2012; Mathews et al. 2012; Casassus 
et al. 2013; Fukagawa et al. 2013; Rosenfeld et al. 2013b; 
Tsukagoshi et al. 2014; Osorio et al. 2014; Huelamo et al. 
2015; Canovas et al. 2015). Some illustrative examples 
are shown in Figure 18 (as well as the insets in Fig. 17). 
Resolved studies of transition disks have so far used hap¬ 
hazard selection, and are severely limited by luminosity 
and resolution constraints. The current sample shows 
cavities with radii of ~15-100AU and emission levels 
suppresed by at least a factor of ~100, surrounded by 
bright, narrow rings (not well-resolved; potentially with 
high optical depths). These morphological features are 
surprisingly more common than the infrared-based tran¬ 
sition disk phenomenon in general, representing a quarter 
or more of the disks in the bright half of the mm luminos¬ 
ity (disk mass) function (Andrews et al. 2011). It is not 
yet clear how much this apparent discrepancy is influ¬ 
enced by the strong selection biases (i.e., bright sources 
only) in the resolved mm continuum sample. 

The simple “dust ring+cavity” model for a transition 
disk structure becomes more complicated when addi¬ 
tional tracers are considered. The presence of a small 
infrared excess (and often silicate emission features) in 
transition disk SEDs demonstrates that the cavities are 
not empty; rather, a “gapped” structure, with an addi¬ 
tional band of small dust grains located near the host 
star, may be a more appropriate description (e.g., Es- 
paillat et al. 2007, 2008, 2010). Since even a low mass of 



Fig. 18.— Examples of the 870 thermal continuum ring mor¬ 
phologies observed for transition disks (at various inclinations) for 
transition disks: ( clockwise, from top left GM Aur (Hughes et al. 
2009b), LkCa 15 (Andrews et al. 2011), RX J1633.9—2442 (Cieza 
et al. 2012), and RX J1604.3—2130 (Mathews et al. 2012). The 
synthesized beams are shown in the lower left of each panel. 

small particles can emit strongly in the infrared, such 
structures can even wash out the characteristic SED 
“dip” signature (e.g., Isella et al. 2010a,b; Andrews et al. 
2011). Moreover, spatially resolved observations of in¬ 
frared scattered light (Dong et al. 2012; Hashimoto et al. 
2012; Follette et al. 2013; Rapson et al. 2015; see Fig. 17 
inset) and sub-mm emission (e.g., Pinilla et al. 2015) find 
smaller particles interior to the mm continuum rings, in 
some cases nearly filling the cavities. 

These trace populations of small particles are clearly 
accompanied by gas. The accretion rates onto the tran¬ 
sition disk hosts tend to be comparable to or slightly be¬ 
low those for more typical disks (Najita et al. 2007, 2015; 
Fang et al. 2009; Espaillat et al. 2012; Ingleby et al. 2014; 
Manara et al. 2014), suggesting that the mass flow is rel¬ 
atively uninterrupted by the gap(s) or cavity. Supple¬ 
mentary signatures of this gas reservoir are found in in¬ 
frared rovibrational lines of common molecules located at 
few AU-scales (e.g., Salyk et al. 2009), unresolved emis¬ 
sion in the wings of millimeter rotational lines (Rosen¬ 
feld et al. 2012), as well as direct spectral imaging of 
those same lines (Casassus et al. 2013; Fukagawa et al. 
2013; Bruderer et al. 2014; Zhang et al. 2014; Perez et al. 
2015; Canovas et al. 2015; van der Marel et al. 2015). 
As with the small grains, constraints on the gas densi¬ 
ties in the cavity are rudimentary. The bulk of the gas 
mass remains confined to the extended structures that 
reach well-beyond the mm continuum rings (Hughes et al. 
2009b; Rosenfeld et al. 2013b; Walsh et al. 2014). 

Considering all these observations together, there is 
strong empirical evidence supporting the theoretical 
model that transition disk structures are largely con¬ 
trolled by a narrow, radial (ring-like) gas pressure max¬ 
imum located outside a zone of substantially depleted 
densities. For the specific examples that have been stud¬ 
ied in detail so far, the measured accretion rates and 
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Fig. 19. — An especially dramatic example of substructure in the 
SAO 206462 transition disk, (left) The 2.2 pm (polarized) scat¬ 
tered light shows prominent spiral structure (Garufi et al. 2013). 
(right) The thermal continuum emission at 435 pm reveals a strong 
azimuthal asymmetry in the ring, which may be related to an un¬ 
resolved spiral pattern (Perez et al. 2014). Both images are on the 
same scale (200 AU on a side), with resolutions on the lower left. 

cavity contents indicate that the pressure gradient out¬ 
side the cavity is not too steep. Such properties are 
generally inconsistent with the predictions of photoevap¬ 
oration models (Alexander & Armitage 2009; Owen & 
Clarke 2012), and so the operating hypothesis is that 
the observed structures are produced by dynamical in¬ 
teractions with very faint (e.g., Pott et al. 2010; Kraus 
et al. 2011), presumably planetary-mass, companions (al¬ 
though more massive stellar companions must be ruled 
out; e.g., Ireland & Kraus 2008). That is certainly an 
exciting prospect, and is naturally driving rapid develop¬ 
ments in the pursuit for details on these specific disks. A 
few putative planetary-mass candidates have even been 
identified inside transition disk cavities (e.g., Huelamo 
et al. 2011; Kraus & Ireland 2012; Reggiani et al. 2014; 
Biller et al. 2014; Quanz et al. 2014). That said, it is 
important to recognize that the sub-population of disks 
receiving most of the attention is still biased (in terms of 
mm luminosity, cavity size, etc.): a larger swath of the 
disk population may have its properties shaped more by 
other evolutionary factors (e.g., photoevaporation). 

6.2.3. Substructure 

A particularly active topic of research on these disks 
is the characterization of their substructure, deviations 
from the simple one-dimensional (radial) description out¬ 
lined above. From a phenomenological perspective, three 
basic types of substructure are observed in transition 
disks: azimuthal asymmetries, spirals, and warps. A 
prominent example (for two of these) is shown in Fig¬ 
ure 19. The former type is common in the limited sam¬ 
ples available, with the millimeter emission rings in some 
cases appearing lopsided: azimuthal contrast levels range 
from modest (a factor of <2; Brown et al. 2009; Isella 
et al. 2013; Rosenfeld et al. 2013b; Perez et al. 2014) 
to extreme (a factor >30; Casassus et al. 2013; van der 
Marel et al. 2013; Fukagawa et al. 2013), and the asym¬ 
metries themselves can account for ^30-100% of the inte¬ 
grated flux. These asymmetric features are azimuthally 
resolved (spanning tens of degrees), but radially nar¬ 
row (<10-20 AU). The currently favored explanation for 
such substructure is the concentration of solids in an az¬ 
imuthal gas pressure enhancement (Birnstiel et al. 2013), 
presumably driven by vortex instabilities produced exte¬ 
rior to the orbit of a perturbing companion (e.g., Barge 
& Sommeria 1995; Klahr & Henning 1997; Wolf & Klahr 


2002; Regaly et al. 2012; Lyra & Lin 2013; Zhu et al. 
2014; Zhu & Stone 2014). 

However, vortices are not the only viable explana¬ 
tions for the observed asymmetries; in some cases, they 
may alternatively be generated by global gravitational 
modes (Mittal & Chiang 2015), eccentricity (Kley & 
Dirksen 2006), or unresolved spiral structures (e.g., Perez 
et al. 2014). The latter is of special interest, since re¬ 
markable spiral patterns have been observed in (usually 
polarized) scattered light from many of the transition 
disks with early-type hosts (a sensitivity-related selection 
bias). Both well-ordered, open spirals (e.g., Muto et al. 
2012; Grady et al. 2013; Garufi et al. 2013) and more 
complex, tightly-wound patterns (e.g., Clampin et al. 
2003; Fukagawa et al. 2004, 2006; Canovas et al. 2013; 
Avenhaus et al. 2014) have been identified, typically ex¬ 
tending from the inner working angles (~10AU) out to 
>100 AU scales. Some debate remains about the origins 
of these patterns, as true density waves (Dipierro et al. 
2014, 2015) or manifestations of vertical substructure 
(Juhasz et al. 2014), but the forthcoming resolution im¬ 
provements at millimeter wavelengths should settle the 
issue. Some care has to be taken in interpreting scat¬ 
tered light images, since the illumination pattern of the 
outer disk surface can also be affected by warped geome¬ 
tries (e.g., Roberge et al. 2005; Quillen 2006; Hashimoto 
et al. 2011; Marino et al. 2015). Resolved spectral line 
data bolsters the evidence that such vertical substruc¬ 
ture may be a common feature in transition disks (e.g., 
Rosenfeld et al. 2012, 2014; Perez et al. 2015). 

With new access to ALMA and high quality (polarized) 
scattered light data, the past few years have seen rapid 
development in resolved studies of transition disks. The 
new signatures of dynamical substructure in these disks 
are marshaling a major investment in both observations 
and theoretical work. Ultimately, the goal is to provide 
firm, predictive links between the observed disk struc¬ 
tures and the properties of the planetary perturbers that 
shape them, with results that have fundamental impacts 
on models of early planetary system architectures, migra¬ 
tion, and the planetary accretion process. At the same 
time, it is beneficial to understand that the “sample” be¬ 
ing studied so extensively is certainly biased: after all, 
these accreting, nun-bright disks with large and partially- 
filled cavities are but a small sub-group of the general 
transition disk population (e.g., Cieza et al. 2008). It 
will be interesting to see how the description of transi¬ 
tion disk structures changes when a more diverse sample 
is analyzed using data of similar quality. 

6.3. Synopsis 

• Disk solids grow and migrate toward the local max¬ 
imum in the gas pressure distribution. The result 
is a vertically and radially size-sorted density struc¬ 
ture, with higher concentrations of larger particles 
at the midplane and in the inner disk. The vertical 
stratification can be observed in the infrared SED, 
or through multiwavelength resolved scattered light 
images. The radial segregation is noted in the spa¬ 
tially resolved mm/radio continuum “colors”. 

• Additional evidence for the inward migration of 
particles due to radial drift is available in the form 
of a steep drop in the dust-to-gas ratio, observa- 
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tionally manifested as a sharp-edged continuum 
profile and a substantial size discrepancy between 
resolved tracers of gas and mm/cm particles. 

• The population of “transition disks” is character¬ 
ized by a ring-like dust geometry, with larger par¬ 
ticles (observed in the mm/cm continuum) concen¬ 
trated in a narrow annulus outside a central, de¬ 
pleted cavity. In the specific examples studied in 
detail so far, gas and smaller particles extend into 
the cavity and out to larger radii. The physical 
interpretation of these structures is related to the 
presence of a radial maximum in the gas pressure 
profile, which effectively traps large particles. 

• A subset of the transition disk population also ex¬ 
hibits a variety of substructures, in the forms of 
azimuthal asymmetries, spiral patterns, and verti¬ 
cal warps. The current thinking is that these mas¬ 
sive, still-accreting systems represent young (giant) 
planetary systems during their formation epoch. 

Additional Reading : the sequence of reviews on the evo¬ 
lution of disk solids by Beckwith et al. (2000), Natta 
et al. (2007), and Testi et al. (2014); the more general 
review on disk evolution by Alexander (2008); and the 
recent overviews of disk dispersal and transition disks by 
Alexander et al. (2014) and Espaillat et al. (2014). 

7. FUTURE DIRECTIONS 

The intentions of this review were to highlight some 
of the fundamental aspects of this field and their key 
applications, and to hopefully serve as a pedagogical re¬ 
source (or starting point) for those interested in joining 
the related research efforts. The past few years have seen 
dramatic advances in the observational study of circum- 
stellar disks, and particularly their connections to the 
formation epoch for planetary systems. It is not a stretch 
to argue that at this stage the subject could be better 
characterized as “observational planet formation”. The 
field is poised for substantial continued development, ow¬ 
ing particular to the capabilities of the nearly completed 
ALMA project. Given that promise, there seems little 
doubt that the landscape of disk research will look sig¬ 
nificantly different in the next decade. Although it is 
difficult to predict that future, the first steps along the 
path toward it are now relatively clear. 

Prominent among these is an effort to sort out some of 
the basic demographics of disk properties, and thereby 
the key factors that impact planet formation efficien¬ 
cies. Section 4 introduced the early efforts along these 
lines, pointing out the preliminary indications that disk 
masses (and thereby the likelihood for planet formation) 
depend on the stellar host mass (Sect. 4.1), multiplicity 
(Sect. 4.2), and environment (Sect. 4.3). But so far the 
limitations on the sizes and properties of samples has sti¬ 
fled further inquiry along these lines. That will change 
with ALMA; several survey projects are underway, and 
more will surely follow. There are many opportunities 
to address some fundamental questions in the field with 
creative explorations of the data available from such sur¬ 
veys. How does the disk mass distribution evolve? What 
are the inter-dependences between and evolutionary be¬ 
havior of the demographic trends that have been noted so 
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Fig. 20.— The HL Tau disk observed at high angular resolution 
(O'/OS, or ~7 AU) in the 870 /im continuum. The series of concentric 
rings are presumably responsible for slowing or halting radial drift, 
and are likely the active, telltale signatures of planetary system 
formation (ALMA Partnership et al. 2015). 

far? What accounts for the tremendous scatter around 
these trends? How does the gas reservoir, as traced by 
key spectral line ratios (e.g., Williams & Best 2014), be¬ 
have with respect to the solids? How are these simple, 
compound diagnostic relations tied to more detailed, el¬ 
emental (i.e., resolved) measurements? 

The last of these questions connects back to the funda¬ 
mental goal of constraining disk density structures (Sec¬ 
tion 5). In the near term, progress toward that goal 
from continuum data alone is a major challenge. Intrin¬ 
sic uncertainties in the opacities of disk solids (Sect. 5.3), 
coupled with the complex processes of growth and mi¬ 
gration that modify their sizes and spatial distributions 
(Sect. 6.1), make a robust and quantitative characteriza¬ 
tion of the densities prohibitive without additional infor¬ 
mation. One way around this impasse is to forge closer 
links to the wealth of spectral line measurements avail¬ 
able from ALMA observations. At this stage, develop¬ 
ing a robust forward-modeling approach to “fit” a large 
and complex union of datasets with an ever-expanding 
(and not well-justified) set of (pseudo-)physical parame¬ 
ters is intractable. Instead, emphasis should be placed on 
identifying and characterizing patterns among resolved, 
multi-tracer diagnostics. For spectral line data, these 
might include the locations of condensation fronts (Qi 
et al. 2011, 2013; Mathews et al. 2013), molecular abun¬ 
dance patterns (e.g., Oberg et al. 2010, 2011; see Dutrey 
et al. 2014), or constraints on the vertical temperature 
distribution (Dartois et al. 2003; Rosenfeld et al. 2013a; 
de Gregorio-Monsalvo et al. 2013). How do these and 
other measurements behave as a function of basic demo¬ 
graphic properties? How are they related to the resolved 
morphologies of the continuum emission? 

An important complementary approach involves build¬ 
ing on the patterns already emerging from resolved con¬ 
tinuum datasets, and associating these with the theo¬ 
retical predictions from models of the evolution of disk 
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solids (Sect. 6.1). Does the tentative correlation between 
continuum luminosities and sizes (Sect. 5.2) hold up for 
larger and different samples, and does it depend on other 
demographic properties? How is the radial variation in 
mm/radio “color” or the location and shape of the outer 
edge of the continuum emission (Sect. 6.1) connected to 
factors like the evolutionary state, stellar host, or the 
structure of the gas reservoir? And most pressing, can 
high angular resolution continuum measurements find 
and characterize evidence for the small-scale substruc¬ 
ture presumed to be responsible for mitigating high ra¬ 
dial drift rates, and thereby preserving large particles in 
the outer disk? The stunning recent ALMA image of the 
canonical HL Tau disk, shown in Figure 20, hints that a 
major conceptual evolution related to this last question 
is underway (ALMA Partnership et al. 2015). 

Substructure observed in these disks will be especially 
important in developing a comprehensive model for the 
formation and early evolution of planetary systems. If 
HL Tau is any indication, observations with ALMA 
may find the general disk population riddled with young 
multi-planet systems carving out concentric, apparently 
resonant, gaps in their natal material. Or perhaps the 
larger-scale asymmetries identified for many transition 
disks (Sect. 6.2) indicate that a less orderly phenomenol¬ 
ogy will be more common. In any case, the underlying 
issue on this topic is finding ways to quantitatively link 
the small-scale features that are observed to the physical 
properties of associated planets and the mechanics of the 
disk-planet interactions. How can multi-tracer measure¬ 
ments of (gas and dust) substructure be optimally com¬ 
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can disk observations be exploited to optimize searches 
for young planets? And, when successful, what can the 
combined measurements of planets and their birth reser¬ 
voirs say about the processes of planetary formation, mi¬ 
gration, and accretion? 

In the larger context of piecing together our cosmic 
origins, these are deep, pressing questions. It is remark¬ 
able to think that the key technological resources needed 
to address these questions are now available, and that 
answers (and new questions) are sure to follow. 
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